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ABSTRACT 

(Distribution  Limitation  Statement  No.  2) 


Techniques  are  described  for  simulating  blast-type  transient  surface  loads  of 
nearly  exponential  pulse  shape  having  characteristic  times  (impulse/peak  pressure) 
r  inging  from  10  to  1000  ysec  (the  quasi-impulsive  range  for  cylindrical  shells 
•oout  1  foot  in  diameter) .  Pressure-time  histories  are  measured  at  various 
positions  around  and  along  cylindrical  models  3.5,  6,  and  12  inches  in  diameter. 

A  basic  set  of  loads  is  obtained  consisting  of  two  limiting  pressure  distributions, 
an  asymmetric  distribution  typical  of  side  exposure  to  a  normally  incident  blast 
wave,  and  a  symmetric  distribution  typical  of  nose-on  exposure.  All  of  the  loads 
are  obtained  using  sheet  explosive  charges  of  various  forms,  from  flat  to  semi- 
cylindrical  to  completely  cylindrical  surrounding  the  model  and  flat  charges 
suspended  at  various  standoffs  in  a  shock  tube.  In  support  of  the  experiments, 
the  self-similar  solutions  for  blast  waves  from  intense  explosions  are  used  to 
calculate  the  range  of  sheet  charges  needed  to  produce  loads  of  interest  and 
to  show  that  the  corresponding  spherical  charges  become  much  too  small  and  much 
too  large  for  practical  application  near  the  extremes  in  load  duration.  Approx¬ 
imate  formulas  are  also  derived  (the  Korobeinikov  theory)  for  the  variation  of 
peak  pressure  with  distance  from  plane,  cylindrical,  and  spherical  charges.  The 
range  of  validity  of  the  formulas  extends  from  high  pressures,  where  the  self¬ 
similar  solutions  are  valid,  to  acoustic  shock  pressures.  Experimental  measure¬ 
ments  from  the  present  program  and  from  compiled  blast  data  show  excellent 
agreement  with  the  theory  over  the  entire  range. 
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SECTION  I 


INTRODUCTION 

This  report  is  concerned  with  economical  methods  for  simulating 
transient  surface  pressure  loads,  similar  to  those  produced  by  blast 
waves,  on  cylindrically  shaped  objects.  The  subject  has  wide  appli¬ 
cation  because  these  loads  are  produced  by  a  variety  of  weapon  effects 
and  cylindrical  shapes  appear  in  a  variety  of  structures  (e.g.,  aircraft, 
missiles,  satellites,  and  reentry  vehicles). 

The  work  is  a  direct  extension  of  the  study  of  response  of  reentry 
vehicle-type  shells  to  blast  loads  under  the  HARTS  program.1  There, 
long  duration  loads  were  produced  by  shock  tubes,  and  a  limited  range 
of  moderately  short  duration  loads,  intermediate  between  shock  tube 
and  impulsi^/e  loads,  was  produced  by  explosive  spheres.  However, 
methods  were  not  available  for  producing  the  entire  range  of  inter¬ 
mediate  duration  loads  with  a  minimum  of  explosive  (less  than  100  pounds). 
This  report  describes  techniques  developed  to  cover  this  intermediate 
range . 

1.  Features  of  Weapon  Effect  Loads 

2 

Surface  loads  from  a  variety  of  weapon  effects  have  three  common 
features  that  make  simulation  with  blast  waves  possible:  the  prrssure 
is  applied  by  a  gas,  it  varies  smoothly  around  the  surface,  and  the 
pulse  shape  (time  history)  consists  of  a  sudden  rise  in  pressure  followed 
by  a  decay  to  zero  pressure.  These  are  the  essential  features  of  the 
diffracted  pressures  experienced  on  the  leading  surface  of  a  vehicle  or 
model  enveloped  by  a  blast  wave.  A  complete  description  of  the  pressure¬ 
time  history  at  every  point  on  the  model  is  quite  involved,  but  to  in¬ 
vestigate  structural  response  it  is  usually  sufficient  to  characterize 
the  load  by  the  product  of  a  pressure  distribution  and  a  pulse  shape, 
i.e.,  the  pressure  distribution  can  be  treated  as  the  same  at  every 
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instant  with  the  pressure  varying  simultaneously  around  the  model  at  a 
single  time  history.  Actual  loads  differ  from  this  idealization;  for 
example,  a  blast  wave  takes  a  finite  time  to  sweep  over  the  model  and, 
also,  the  low  pressure  on  the  far  side  of  the  model  rises  slowly  rather 
than  suddenly.  However,  these  complicating  features  are  generally  secon¬ 
dary,  either  because  time  differences  are  small  compared  to  structural 
response  times,  or  because  the  area  in  which  the  pressure  is  not  proper¬ 
ly  represented  by  a  product  form  is  at  a  relatively  low  pressure.1 


The  pressure  distribution  for  a  particular  weapon  effect  depends  on 

the  orientation  of  the  model  relative  to  the  on-coming  load  source.  For 

example,  radiation  induced  surface  loads  on  a  cylinder  illuminated  from 

one  side  produces  a  cos  0  distribution  over  the  illuminated  side,  and 

no  pressure  on  the  shaded  side.  Side  exposure  to  intense  blast  waves 

2 

produces  an  approximately  cos  Q  distribution  over  the  exposed  side 
superimposed  on  a  relatively  low  pressure  that  can  be  treated  as  uniform 
around  the  circumference.  End-on  exposure  results  in  pressures  sym¬ 
metric  around  the  circumference,  while  intermediate  incidence  angles 
result  in  varying  degrees  of  asymmetry.  Two  representative  distribu¬ 
tions  are  sketched  in  Fig.  1.  There  is  no  end  to  the  number  of  distri¬ 
butions  that  can  be  imagined,  but  it  is  reasonable  to  expect  that  the 
distinguishing  feature  of  a  load  distribution  is  its  degree  of  symmetry, 
and  it  is  common  practice  to  divide  weapon  encounters  into  two  basic 
types,  side-on  (asymmetric)  and  nose-on  (symmetric).  Asymmetric  loads 
are  the  most  probable  and  generally  result  in  more  severe  structural 
damage  and  have  therefore  received  the  most  attention. 

Separation  of  surface  loads  into  asymmetric  and  symmetric  types 
is  a  logical  division  for  structural  response.  The  primary  response  of 
a  cylinder  or  cone  to  symmetric  loads  is  the  development  of  hoop  mem¬ 
brane  stresses  and  possibly  buckling.  For  asymmetric  loads,  a  net  side 
load  must  also  be  sustained;  this  gives  rise  to  beam  bending  and  shear 
forces  and  also  to  large  end  resction  forces  at  points  of  attachment, 
possibly  causing  cracks  or  local  yielding.  Accelerations  of  internal 
components  are  also  generally  more  severe  for  asymmetric  loads.  These 
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FiG.  1  REPRESENTATIVE  PRESSURE  DISTRIBUTIONS  AND  PULSE  SHAPE 

observations  suggest  that  the  extremes  of  structural  response  can  be 
examined  by  testing  with  normally  incident  side~on  exposure  and  with 
nose-on  exposure.  Tests  on  cylindrical  shells1  demonstrate  that  the 
exact  pressure  distribution  of  a  side-on  load  is  secondary  and,  in  fact, 
for  shell  buckling,  critical  pulses  for  damage  from  asymmetric  loads  can 
be  nearly  the  same  as  from  symmetric  loads. 


The  essential  features  of  the  pulse  shape  are  exhibited  by  the 
exponential  pulse 

0  t  <  0 

-t/t 

P  e  °  t  >  0 

where  P  is  the  peak  pressure,  t  is  time,  and  t  is  a  characteristic 

time  related  to  the  total  impulse  under  the  pressure-time  curve  by 

t  =  I/P.  A  sketch  of  this  pulse  is  shown  in  Fig.  1.  The  sensible  dura- 
o 

tion  of  the  pulse  is  about  5tQ.  Actual  pulses  differ  from  this  exponen¬ 
tial  shape,  but  the  value  of  5  is  a  representative  ratio  between 
sensible  duration  and  an  artificial  characteristic  time  defined  by  I/P. 
For  either  the  idealized  exponential  or  the  actu&x  pulse,  selection  of 
a  sensible  duration  is  quite  subjective;  to  avoid  this  difficulty  we 
shall  use  instead  the  artificial  characteristic  time  t  =  I/P  for  any 
pulse.  If  the  pulses  do  not  differ  widely  from  a  single  pulse  shape, 
they  can  be  characterized  by  two  parameters,  for  example  peak  pressure 
and  impulse,  or  peak  pressure  and  characteristic  time.  Surface  loads 
from  weapon  effects  of  inters  it  here  have  pressures  ranging  from  many 
kilobars  to  a  few  bars,  and  characteristic  times  ranging  from  a  fraction 
of  a  microsecond  to  many  milliseconds. 

2.  Basis  for  Simulation  and  Structural  Response  Investigations 

To  determine  structural  response  to  each  weapon  effect  separately 
would  be  both  expensive  and  time-consuming.  If  the  problem  is  inverted, 
however,  and  structural  response  is  first  determined  for  a  range  of 
loads  significant  to  the  structure,  it  will  then  be  relatively  easy  to 
determine  structural  response  for  a  particular  weapon  effect.  To  fully 
exploit  this  approach , it  is  essential  that  structural  response  be  under¬ 
stood  over  the  entire  range  of  loads;  only  then  can  the  influence  of 
load  characteristics  on  structural  response  be  discussed  in  terms  of  the 
structure  rather  than  in  terms  of  a  specific  weapon  effect. 
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From  the  preceding  discussion,  and  arguments  given  in  Reference  1, 
we  take  the  essential  surface  load  parameters  to  be  peak  pressure, 
impulse,  and,  to  a  lesser  extent,  the  symmetry  of  the  pressure  distri¬ 
bution.  An  adequate  description  of  the  strength  of  a  given  structure 
for  weapon  effects  work  can  then  be  obtained  by  tests  using  two  pressure 
distributions,  one  symmetric  and  the  other  asymmetric,  and  pressure 
pulses  ranging  in  peak  pressure  and  impulse  over  the  range  of  signifi¬ 
cance  to  the  structure.  All  of  the  remaining  details  of  the  actual 
loads  are  expected  to  be  of  secondary  importance.  These  include  effects 

of  pulse  shape  (e.g.,  triangular  vs.  exponential  decay),  detailed  pres- 

2 

sure  distribution  (e.g.,  cos  9  vs.  cos  0  load  on  one  side,  or  the 
effect  of  small  pressures  on  the  unloaded  side)  and  time  phasing  (e.g., 
inexact  simulation  of  the  time  for  a  blast  wave  to  sweep  over  the  model). 
Since  this  simplification  is  merely  a  reasonable  assumption  and  not  an 
a  priori  conclusion,  the  approach  is  t.o  examine  the  effect  on  structural 
response  of  changing  these  details  within  the  capabilities  of  simple 
tests.  If  the  resulting  effects  are  indeed  small,  then  the  main  bulk 
of  experiments  can  be  run  using  more  standard  tests  with  fewer  parameters 
This  results  in  a  systematic  method  for  determining  structural  vulner¬ 
ability  and  avoids  the  excessive  expense  of  attempting  "exact"  simulation 
which  can  never  be  attained. 

Peak  pressure  and  impulse  are  chosen  as  the  two  parameters  to 
characterize  the  pulse  because  these  are  the  significant  parameters  at 
each  extreme  of  pulse  duration;  for  very  short  durations  response  depends 
only  on  impulse  and  for  very  long  du.'  ations  response  depends  only  on  peak 
pressure.  For  intermediate  values  the  response  depends  on  both  pressure 
and  impulse.  Critical  curves  for  a  given  extent  of  damage  often  have 

r 

a  hyperbolic  shape  in  the  pressure-impulse  plane;  for  example.  Fig.  2 
(taken  from  Reference  1)  shows  curves  for  simple  cylindrical  shells  of 
aluminum  subjected  to  asymmetric  loads.  In  this  plane,  lines  of  constant 
characteristic  time  t  are  rays  through  the  origin  or,  in  the  log-log 
plot  of  Fig.  2,  parallel  lines  at  unit  slope.  To  experimentally  deter¬ 
mine  critical  damage  curves,  methods  must  be  available  to  provide  pulses 
throughout  the  entire  region  of  interest  in  the  pressure-impulse  plane 
for  both  asymmetric  and  symmetric  loads. 
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FIG.  2  PRESSURE-IMPULSE  DIAGRAM  FOR  SEVERAL  SIMULA TIGN  METHODS 

(for  6-inch-diameter  cylinders;  hyperbolic  curves  ore  buckling  thresholds 
for  6061“ T6  aluminum  shells,  L/D  =  1) 


3.  Methods  of  Surface  Load  Simulation 

The  range  of  application  of  the  simulation  methods  used  in  the 
previous  programs  is  also  shown  in  Fig.  2,  for  6-inch-diameter  models. 
Loads  in  the  region  of  very  short  pulses,  to  the  left  of  the  t  = 

O-Oi  nsec  line,  are  obtained  with  explosive  sheet  or  gas  detonated 
adjacent  to  the  model.  *  *  For  loads  of  intermediate  duration  the 
two  steep  curves  on  the  left  give  the  peak  pressure  and  impulse  at  the 
point  of  maximum  pressure  for  side  exposure  to  blast  waves  from  0.5  and 
4-pound  charges.  The  parameter  varied  on  each  curve  is  the  distance 
(standoff)  from  the  model  to  the  charge.  For  long  duration  loads  the 
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three  curves  at  approximately  45  on  the  right  are  for  side  exposure  to 

X  S  6 

blast  waves  from  explosively  driven  shock  tubes  ’  ’  10,  30,  and  100 

feet  iong.  The  parameter  varied  on  these  curves  is  the  weight  of  ex¬ 
plosive  in  the  shock  tube. 

The  object  of  the  present  investigation  is  to  provide  better  load¬ 
ing  methods  in  the  intermediate  region  between  contact  loads  and  shock 
tube  loads.  Characteristic  times  in  this  region  range  over  0.01  <  t  < 

1  msec,  a  factor  of  100.  Since  the  characteristic  time  for  a  given 

pressure  varies  in  direct  proportion  to  the  charge  radius  and  hence  as  the 

* 

cube  root  of  spherical  charge  weight,  spherical  charges  ranging  in 

g 

weight  by  a  ratio  of  more  than  10  would  be  required  to  cover  this  region. 
For  example,  at  surface  pressures  near  2000  psi,  from  Fig.  2  a  one-pound 
sphere  gives  t  =  0.03  msec.  To  obtain  t  -  1  msec  would  require  a 
charge  weighing  more  than  (1/0. 03)3  =  37,000  pounds,  much  too  large  for 
economical  use.  At  the  other  extreme,  to  obtain  a  pressure  pulse  at 
1000  psi  and  the  lower  limit  of  interest,  t  =  0.01  msec,  would  require 
a  charge  of  0.01  pound,  much  too  small  for  structural  testing.  Further¬ 
more,  the  rapid  variation  of  peak  pressure  with  standoff  from  a  spherical 
charge  makes  accurate  adjustment  of  pressure  difficult. 

All  of  these  shortcomings — excessively  large  charges  and  charges  so 
small  that  they  do  not  provide  coverage  of  the  model,  and  extreme  accuracy 
requirements  on  standoff  distance — can  be  overcome  by  using  sheet  charges. 
In  their  simplest  form  these  consist  of  a  flat  sheet  of  explosive  suspended 
in  air  as  for  spherical  charges.  Unlike  the  sphere,  however,  at  distances 
from  the  sheet  small  compared  to  its  lateral  dimensions  the  blast  wave 
expands  in  only  one  dimension  rather  than  in  three,  '"bus  for  a  fixed  area, 
charge  weight  varies  with  the  first  power  of  the  characteristic  time  rather 


This  simple  geometric  scaling  applies  to  the  pressure  pulse  incident 
upon  the  model.  For  very  large  charges  the  impulse  experienced  on  the 
model  would  be  smaller  than  given  by  this  rule  because  of  flow  around 
the  model,  thus  requiring  a  still  larger  range  of  charge  weights. 

The  term  sheet  charge  is  used  to  designate  any  charge  in  the  form  of  a 
two-dimensional  surface.,  whether  it  is  flat,  curved,  or  in  the  form  of 
a  complete  cylinder.  It  may  or  may  not  be  made  of  sheet  explosive,  a 
name  reserved  here  for  plasticized  PETN  explosive,  made  by  Du  Pont  under 
the  registered  trade  names  Detashcet  C  ar.d  Detasheet  D. 
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than  the  third,  thereby  reducing:  the  required  weight  ratio  from  10  to 
2 

10  .  Blast  loads  of  interest  can  be  obtained  with  average  explosive 
sheet  thicknesses  ranging  from  a  few  mils  to  something  less  than  an 
inch.  To  maintain  the  lateral  dimensions  of  the  charge  large  compared 
to  the  standoff  requires  excessively  large  charge  areas  for  distances 
greater  than  a  few  feet.  At  these  larger  distances  the  charge  can  be 
placed  in  a  shock  tube  to  maintain  cne-dimensional  flow.  At  very  small 
distances  the  charge  must  be  curved  partly  around  the  model  to  obtain 
an  appropriate  pressure  distribution  and  loading  simultaneity;  in  the 
limit  it  becomes  a  contact  charge.  Thus  sheet  charges  continuously  span 
the  load  region  of  interest  and  at  each  extreme  merge  continuously  with 
contact  charge  and  shock  tube  techniques. 

A  pressure-impulse  diagram  similar  to  Fig.  2  can  also  be  drawn  for 
symmetric  loads.  In  this  case  the  sheet  charges  (both  contact  and 
standoff)  take  the  form  of  cylinders  completely  surrounding  the  model, 
and  in  the  shock  tube  the  model  is  pointed  nose-on  down  the  tube  (with 
a  suitable  nose-tip,  if  not  already  part  of  the  model).  For  either  sym¬ 
metric  or  asymmetric  loads  on  models  about  one  foot  in  diameter,  the  char¬ 
acteristic  times  which  (somewhat  arbitrarily)  divide  contact  charges,  sheet 
blast  charges,  and  shock  tube  techniques  are  20  Usee  and  1  millisecond.  A 
schematic  diagram  of  the  corresponding  load  regions  is  shown  in  Fig.  3. 
Contact  techniques  are  limited  on  the  left  by  minimum  practical  sheet 
or  gaseous  explosive  thicknesses,  and  shock  tube  techniques  are  limited 
on  the  right  by  the  economics  of  making  very  long  shock  tubes. 
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SECTION  II 


SUMMARY 


1.  Theory 

Using;  the  self-similar  theory  of  intense  explosionss  it  is  shown 
in  Section  III  that  blast  waves  from  flat  sheet  charges  have  free-field 
pressure-time  histories  very  similar  to  those  from  spherical  charges. 
Thus,  at  these  high  overpressures,  sheet  charges  give  a  good  simulation 
of  spherical  blast  waves.  The  theory  is  also  used  to  make  a  quantita¬ 
tive  congjarison  between  the  explosive  required  in  the  two  geometries. 

At  lower  pressures,  further  from  the  explosion,  the  self-similar 
theory  underestimates  the  pressure  because  the  energy  of  the  ambient  gas 
enveloped  by  the  blast  wave  is  neglected.  Using  an  approximation  sug¬ 
gested  by  Korobeinikov,  it  is  shown  that  the  decay  of  peak  pressure  with 
distance  from  the  charge,  taking  into  account  the  initial  pressure  p^ 
of  the  ambient  gas,  can  be  expressed  by  very  simple  formulas  for  all 
three  geometries.  Peak  incident  shock  overpressure  is  given  by 


-  _si_A2r1  + <x+a-v/2] 

P,  Y  +  1  L  J 


where  y  is  the  specific  heat  ratio  of  the  ambient  (ideal)  gas  and  A 
is  given  by  the  following  formulas  in  each  geometry: 


Plane  Waves 

A  - 

9<Ti* 

all  B 

2 

-2 

i6ayR 

R  £  2 

Cylindrical  Waves 

A  - 

16  V2"  GtYR3/2 

R  £  2 

- 

-2 

j  25cyR3 

R  £  2 

Spherical  Waves 

A  = 

50avR2  +  3] 

R  2:  2 

(2) 
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In  each  case  R  is  a  dimensionless  radius  defined  in  terms  of  the  actual 
radius  (standoff)  x  by 


where  V  =  1,  2,  or  3  in  plane,  cylindrical,  or  spherical  geometry,  re¬ 
spectively,  and  Eq  is  the  energy  deposited  in  the  blast  wave  by  the 
expletive,  given  approximately  by  the  yield  energy  of  the  explosive  per 
rnit  area,  per  unit  length,  or  total  energy,  for  v  =  1,  2,  or  3.  The 
parameter  a  in  Eq.  (2)  is  a  function  of  y  in  each  geometry;  for  air, 
y  =  1.4  and  a  =  1.075,  1.C0,  and  0.85  for  v  =  1,  2,  3.  Taking 

E  =  1430  cal/gm  for  50/50  oentolite,  Eq.  (2)  gives  pressures  within 
o 

5  percent  of  compiled  experimental  data  for  spherical  charges  to  pressures 

as  low  as  (p  /p,  )  -  1  =  0.03.  Good  agreement  is  also  fou.*d  for  the 

plane  blast  waves  measured  here,  but  measured  pressures  extend  only  as 

low  as  (p„/p  )  -  1  «  2. 

^  1 

2.  Simulation  of  Asymmetric  Loads 

The  desired  range  of  quasi-impulsive  loads  was  obtained  with 
(1)  flat  sheet  charges  placed  on  the  ground,  with  the  model  suspended 
above  them,  and  (2)  flat  sheet  charges  placed  in  a  2-foot-diameter  shock 
tube  with  the  model  at  the  open  end  of  the  tube.  Standoff  distances 
were  varied  from  a  few  inches  to  20  feet,  and  the  sheet  charge  thick¬ 
ness  was  varied  from  a  few  mils  to  a  few  tenths  of  an  inch.  Pressure- 
impulse  curves  for  these  techniques  are  given  in  Fig.  4.  The  three 
uppermost  curves  are  for  charges  on  the  ground  at  standoffs  of  0.5,  1, 
and  2  feet,  and  the  four  lower  curves  are  for  charges  in  the  shock  tube 
at  standoffs  of  2.5,  5,  10,  and  20  feet.  The  parameter  varied  along 
each  curve  is  the  density  (or  equivalent  thickness)  of  the  charge  in 
pounds  per  square  foot  and  in  mils. 


* 

A  mixture  of  PETN  and  TNT. 


11 


PEAK  PRESSURE 


FIG.  4  PRESSURE-IMPULSE  CURVES  FOR  ASYMMETRIC  LOADING  TECHNIQUES 
(6-inch-diom«t«r  modal) 


Loads  were  also  determined  for  larger  standoffs  from  charges  on  the 
ground,  but  the  same  loads  can  be  obtained  with  less  explosive  by  using  the 
shock  tube.  Similarly,  standoffs  somewhat  smaller  than  2.5  feet  can  be 
used  in  the  shock  tube,  but  greater  complications  arise  in  initiating 
the  charge  at  many  points  to  obtain  a  nearly  plane  blast  wave.  (This 
difficulty  is  avoided  for  charges  on  the  ground  by  using  a  running 
detonation,  which  is  acceptable  because  there  are  no  tube  walls  from 
which  the  resulting  oblique  shock  can  reflect  to  give  undesirable  second¬ 
ary  shocks  on  the  model.)  The  crossover  standoff  of  about  2  feet  between 
the  two  techniques  is  satisfactorily  within  these  limitations  for  the 
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i-inch-diameter  model  and  2-foot-diameter  shock  tube  used  here.  Impulses 
at  the  same  sttndoff  are  about  40  percent  higher  in  the  shock  tube  be¬ 
cause  there  is  no  loss  from  lateral  expansion. 

Comparison  of  the  measured  peak  pressures  with  the  self-similar 
theory  shows  that  over  the  entire  range  of  loads  the  peak  pressure  is 
the  reflected  pressure  that  would  be  experienced  on  a  flat,  rigid  wall. 
Pulse  shape  is  approximately  exponential  as  in  Fig.  1.  At  the  larger 
standoffs,  however,  an  exponential  fit  to  the  main  body  of  the  pulse 
gives  a  peak  pressure  of  only  about  70  percent  of  the  measured  initial 
pressure,  possibly  because  the  duration  of  diffraction  flow  is  short 
compared  to  the  total,  pulse  duration.  The  distribution  of  peak  pressure 
around  the  front  half  of  the  model,  facing  the  blast,  was  found  to  be 
given  within  ±10  percent  by 


P 


<P  - 
r 


Pi>  cos  8  + 


90°  <  Q  <  90°  (4) 


where  P^  and  P^  are  the-  reflected  and  incident  overpressures.  Few 
measurements  were  made  on  the  back  half  of  the  model  because  the  pres¬ 
sures  were  very  low  compared  to  the  front  face  pressure  and  mechanical 
gage  vibrations  made  accurate  measurement  difficult.  These  pressures 

were  always  less  than  P.  and  ,  since  P  >  >  P.  for  loads  of  interest, 

l  r  i 

were  considered  to  be  of  little  importance  for  structural  response. 

3.  Simulation  of  Symmetric  Loads 

Figure  5  gives  pressure-impulse  curves  for  symmetric  loads  obtained 
using  (1)  unconfined  cylindrical  charges  and  (2)  flat  charges  in  the 
shock  tube  with  the  model  facing  axially  down  the  tube  and  backed  by  a 
plate  so  that  the  model  experiences  the  reflected  pressure  propagating 
from  the  back  plate.  The  data  are  less  extensive  than  for  asymmetric 
loads  because  these  techniques  were  previously  untried  and  the  experi¬ 
ment!  were  mainly  exploratory.  The  upper  two  curves  show  that  cylindri¬ 
cal  charges  2  feet  and  4  feet  in  diameter  give  loads  having  characteris¬ 
tic  times  near  t  =50  and  100  Psec.  The  lowermost  curve  shows  that 
o 

the  shock  tube  gives  loads  having  tQ  an  1500  Psec  when  the  charge  is 
placed  at  the  base  of  the  tube  (standoff  x  =  18,5  feet). 
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FIG.  5  PRESSURE-IMPULSE  CURVES  FOR  SYMMETRIC  LOADING  TECHNIQUES 

(6-inch-diam«t*r  model) 

Satisfactory  loads  having  t  w  500  M-sec  could  probably  be  obtained 
by  placing  the  charge  at  x  =  5  feet  in  the  shock  tube,  but  still  smaller 
standoffs  could  give  objectionable  shock  decay  along  the  length  of  the 
model  and  durations  of  the  incident  pressure  pulse  comparable  to  that 
of  the  reflected  pulse  from  the  model  back  plate.  This  lower  limit  of 
pulses  obtainable  in  the  capped  shock  tube  technique  is  indicated  by  the 
dashed  curve  in  Fig.  5.  Pulse  times  from  t  -  100  to  500  Psec  would 
require  cylindrical  charge  diameters  larger  than  4  feet,  which  may  give 
difficulty  in  maintaining  the  stability  of  the  imploding  shock.  Even 
at  the  4-foot-diameter,  sigr.ificant  variations  in  peak  pressure  around 
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the  model  were  observed  in  the  highest  charge  density  shot.  A  continuing 
experimental  and  theoretical  study  would  be  required  to  give  a  range  of 
symmetric  loads  comparable  to  the  range  of  asymmetric  loads  demonstrated 
in  Fig.  4. 

4.  Variations  in  Load  Distribution  and  Pulse  Shape 

A  series  of  experiments  were  also  performed  to  explore  possible 
methods  for  varying  load  distribution  and  pulse  shape.  To  vary  load 
distribution,  the  technique  of  varying  charge  density  around  the  model, 
which  had  already  been,  demonstrated  to  provide  variations  in  impulse 
distribution  for  contact  charges,  was  applied  to  sheet  charges  at  stand¬ 
offs  of  a  few  feet.  It  was  found  that  charge  density  variations  in  a 
flat  charge  had  little  effect  on  pressure  distribution.  Density  varia¬ 
tions  in  a  curved  charge,  however,  had  a  strong  effect  on  pressure  dis¬ 
tribution.  In  fact,  changes  in  curvature  (from  flat  to  semi-cylindrical 
around  the  model)  at  a  constant  charge  density  also  gave  significant 

changes  in  pressure  distribution.  Using  both  flat  and  curved  charges 

2 

pressure  distributions  from  cos  0  to  uniform  over  the  loaded  side  were 
obtained. 

Changes  in  pulse  shape  were  studied  for  shock  tube  loads  using 
volume  charges  of  various  lengths  and  density  distributions.  Flat 
charges  (as  used  for  Fig.  4)  were  found  to  give  the  most  satisfactory 
simulation  of  expected  high  altitude  pressure  pulses.  Charges  having  a 
stepwise  variable  density...  most  dense  near  the  model,  gave  pulse  shapes 
approaching  that  of  the  flat  charge  but  had  undesirable  pressure  jumps, 
probably  caused  by  the  discontinuities  in  charge  density.  Such  charges 
should  be  used  only  if  the  stresses  in  the  shock  tube  from  the  localized, 
more  intense  flat-charge  explosions  would  cause  tube  damage. 


Section  III 


SIMULATION  OF  ASYMMETRIC  LOADS 


Asymmetric  loads  are  more  probable  and  tend  to  be  more  damaging 
than  symmetric  loads  and  are  therefore  treated  first  and  more  extensively 
The  theoretical  feasibility  of  using  sheet  chsyrges  is  demonstrated  and 
then  extensive  measurements  .of  pressure,  impulse,  and  pressure  distribu¬ 
tion  around  cylindrical  models  are  given  for  free  expansion  from  flat 
charges  and  for  flat  charges  suspended  in  a  twc-foot-diameter  shock 
tube. 


1.  Theoretical  Pressure  Pulses 

Before  describing  the  experimental  results  it  is  helpful  to  make  a 
theoretical  estimate  of  sheet  charges  required  for  the  desired  pressure 
pulses  and  compare  them  with  the  corresponding  spherical  charges.  This 
can  be  done  using  the  theory  of  blast  waves  from  intense  explosions 
given  :ln  Section  VI,  Near  the  explosion  the  self-similar  theory  gives 
the  following  expressions  for  peak  incident  pressure  pi  and  time  scale 


8R 


-v 
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(6) 


where  p^  and  are  the  pressure  and  density  of  the  undisturbed 

ambient  gas  ahead  of  the  shock,  y  is  the  specific  heat  ratio,  v  is 
an  index  having  values  v  =  1,  2,  and  3  in  plane,  cylindrical,  and 
spherical  blast  waves,  and  is  the  corresponding  constant.  The 

distance  from  the  center  of  the  charge  to  the  shock  front  is  denoted 
here  by  x  and  in  Eqs.  (5)  and  (6)  appears  in  the  dimensionless  form 


R  ,  normalized  to  a  characteristic  length  in  each  geometry  according  to 
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In  each  geometry  E^  is  the  effective  release  energy  of  the  charge, 
having  units  of  energy  per  unit  area  in  plane  geometry,  energy  per  unit 
length  in  cylindrical  geometry,  and  total  energy  of  the  charge  in  spheri¬ 
cal  geometry. 

a.  Comparison  of  Spherical  and  Plane  Blast  Waves 

We  first  compare  the  incident  pressure  pulse  shapes  from  plane 
and  spherical  charges.  These  are  shown  in  Fig.  6,  taken  from  Fig.  52a 
with  the  time  scale  for  the  plane  wave  compressed  by  a  factor  of  2  to 
give  close  agreement  at  early  times.  With  this  modification, the  two 
pulses  are  nearly  identical  until  p/p±  falls  below  0.2.  At  later 
times  the  pressure  decays  more  slowly  for  the  plane  wave,  but  the  dif¬ 
ference  is  not  great  and  effects  not  treated  in  the  self-similar  theory 
would  tend  to  make  the  pulses  more  nearly  the  same.  Similar  comparisons 
of  the  density,  particle  velocity,  and  dynamic  pressure  behind  the  shock 
can  also  be  made;  these  again  show  some  differences  between  the  two 
waves  but  the  similarity  is  close  enough  that  plane  blast  waves  give  a 
good  simulation  of  spherical  blast  waves. 

Knowing  that  the  pulse  shapes  are  similar,  we  can  now  compare 
the  charge  weights  and  standoff  distances  for  plane  and  spherical 

blast  waves  that  give  the  same  peak  pressure  and  duration.  Equating 
peak  pressures  from  Eq.  (5)  gives 

25a3  3  3 

R1  =  -^*3=2-16R3-  (8) 

in  which  =  1.075  and  =  0.85  have  been  taken  for  y  =  1.4  as 


The  notation  here  is  slightly  changed  from  that  in  Section  VI  because 

wayes  of  different  geometries  are  being  compared:  a  subscript  v  is 

used  here  and  r  is  replaced  by  s  to  avoid  confusion, 
o  v 
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t/ff  {SPHERE} 
t/2!f  (PLANE) 

FIG.  6  INCIDENT  PRESSURE  PROFILES  FROM  SELF-SIMILAR  THEORY 


discussed  in  Section  VI.  Fro*  Fig.  6  the  durations  are  made  comparable 
by  taking  2t{1  =  tf3<  Using  this  in  Eq.  (6)  results  in 


1/2  _  _  3/2  1/2  „  5/2 


*°T  S1 


«3  S3  *3 


Combining  this  with  Eq.  (8)  yields 


s 

1  2  \125  ) 


is  R  -2  _  0.137  3 
O  3  -~2“*3 

^  X3 


Frca  the  def.'nition  of  s^  in  Eq,  (7),  s1  is  proportional  to  the 
thickness  of  the  sheet  charge  (weight  per  unit  area)  and  s_  is  pro- 

-  O  ^ 

portions.!  to  the  radius  of  the  spherical  charge  (s„  is  proportional  to 


the  total  charge  weight),  in  these  terns  Eq.  (10)  becomes 


0.137  w 

W1  =  T“  *3’ 


or 


W3  =  81  Vl 


(11) 


where  w  is  weight  per  unit  area  of  the  sheet  charge  and  W  is  total 
weight  of  the  spherical  charge.  Conbining  Eqs.  (8)  and  (10)  with  the 
definition  of  from  Eq.  (7)  gives 


x 


1 


10  X3 


(12) 


Equations  (11)  and  (12)  are  the  desired  relations  between 
charge  sizes  and  model-to-charge  distances  and  are  applicable  for 
pressure  pulses  in  the  range  of  approximate  validity  of  the  self-similar 
approximation.  This  range  is  rather  narrow  (5  <  p^/p^  -  1  <  50),  but 
the  corresponding  range  in  reflected  pressure  with  Pj  =  14.7  psi  is 
300  <  pf  <  5000  psi,  which  happens  to  be  the  range  in  which  many  struc¬ 
tural  models  of  interest  are  damaged .  Also,  since  Eqs.  (11)  and  (12) 
are  comparative,  their  range  of  validity  is  probably  somewhat  larger 
than  the  range  of  absolute  accuracy  of  the  self— similar  theory.  From 
Eq.  (12)  are  see  directly  that  the  space  required  for  a  sheet  charge 
experiment  is  inherently  less  than  that  required  for  a  spherical  charge, 
the  standoff  distance  being  only  30  percent  of  that  for  a  spherical 
charge.  From  Eq.  (11)  we  see,  as  expected,  that  at  large  distances 
(long  durations)  a  great  savings  in  weight  can  be  achieved  by  using  one¬ 
dimensional  blast  waves. 

b.  Estimates  of  Sheet  Thickness  and  Standoffs 

As  discussed  in  Section  I,  the  model  pulse  characteristic  times 

for  which  sheet  charges  appear  attractive  are  in  the  range  10-5  <  t  < 

-3  ° 

10  sec.  Using  the  self-similar  theory,  the  range  of  sheet  charge 

thicknesses  and  standoff  distances  required  to  give  these  characteristic 

times  for  the  incident  pressure  pulses  can  be  estimated.  The  theory  will 

also  prove  valuable  later  In  interpreting  peak  pressure  measurements  in 

the  experiments. 


tioas  csh 


where  h  is  the  explosive  thickness  and  Y  is  its  effective  release 

energy  per  unit  volume.  PETN  has  an  energy  release  of  about  1600  cal/gs 

3  11  3 

which,  at  a  1.4  gm/cm  density,  gives  Y  a#  10  ergs/cm  .  Substituting 

%  - 

Eq.  (13)  into  Eq.  v5)  gives  for  the  explosive  thickness  ,  • 

kl/2 


!i  (1 lV 

*  Wi/ 


3/2  tfl 


(14) 


CoHparison  of  Figs.  1  and  6  shcvs  that  to  a  reasonable  approximation 

t  for  plane  blast  waves  can  be  used  as  the  equivalent  exponential 

characteristic  tine  t  .  Thus  we  can  substitute  10  <  t  <  10 

o  o 

directly  into  Eq.  (14)  with  p  =  10  dyne/ca  ,  p  =  1.293  x  10 

3  11  3  1  ^ 

ga/ca  ,  Y  =  10  ergs/cm  ,  and  =  0.012  to  obtain  the  range  for 

h  at  this  pressure: 


0.0020  <  h  <  0.20  cm 


(15) 


To  obtain  the  corresponding  standoff  distances,  we  use  Eqs.  (7)  and  (13) 
to  get 


a*.  =  R  h  ~ 

1  P, 


which,  with  Eq.  (15),  yields 


2.4  <  x  <  240  ca 


(16) 


Since  the  characteristic  tines  of  the  diffracted  pressures 
measured  on  cylindrical  aodels  are  scaewhat  less  than  the  incident 
tines  (about  half,  aud  decreasing  with  increasing  standoff,  as  we  shall 
see)  Eqs.  (15)  and  (16)  give  lower  bounds  on  the  required  values. 

Doubling  Miese  values  we  obtain  standoff  distances  from  about  2  inches 
to  20  feet  and  explosive  thicknesses  fron  about  2  mils  to  0.2  inch. 

These  ranges  can  be  reasonably  obtained  experimentally.  The  lower  limit 
of  2  mils  for  explosive  thickness  cannot  be  obtained  in  a  continuous 
sheet,  but  thin  extruded  rods  of  explosive  can  be  used  to  give  an  average 


thickness  of  2  mils  at  a  spacing  small  compared  to  the  2-inch  standoff. 

At  pressures  appreciably  belovr  the  3000  psi  example  here,  however, 

gaseous  explosives  would  have  to  be  used  to  maintain  durations  as  short 
2  3 

as  10  Usee.  *  The  upper  limit  of  x  =  20  ft  requires  a  shock  tube  to 
maintain  one-dimensional  flow,  as  already  observed  in  Section  I. 

I 

By  comparison,  the  spherical  charges  corresponding  to  the 
range  of  flat  charges  in.Eq.  (15)  and  (16)  become  impractical  at  each 

3 

extreme.  Using  Eqs.  (11)  and  (12)  and  an  explosive  density  of  1.4  gm/cm  , 
the  range  of  spherical  charge  weights  and  standoffs  are 

1.30  <  W  <  1.30  x  106  gm 

«3 

(17) 

8  <  x3  <  800  cm 

2.  Experimental  Arrangements 

a.  Model  Support  and  Charge  Layout 

The  general  arrangement  typical  of  most  of  the  field  experi¬ 
ments  is  shown  in  Fig.  8.  The  charge?  are  layed  on  a  one-inch-thick 
steel  plate  resting  on  a  bed  of  sand.  Considering  the  plate  as  a  rigid 
wall,  this  doubles  the  effective  yield  of  the  explosive  as  compared  to 
expansion  from  a  charge  freely  suspended  in  air-.  The  model  is  held  in. 
a  stand  large  enough  to  be  used  over  charges  12  feet  square.  Standoff 
is  adjusted  by  raising  or  lowering  the  entire  cross  frame.  Coarse  ad¬ 
justment  is  made  by  sliding  the  frame  over  the  four  pipe  legs:  fine 

adjustment  is  made  by  telescoping  screw  jacks  at  the  base  of  each  leg. 

* 

The  charges  in  Fig.  8  are  all  made  of  Primacord,  which  consists 
of  PETH  packed  in  a  polyethylene  sheath  and  textile  braid.  Primacord 
was  used  for  most  of  the  shots  because  it  is  relatively  inexpensive 
(about  $3  per  pound  of  core  loading)  and  is  readily  available  from  local 
suppliers  in  sizes  from  25  to  400  grains  of  PETN  per  foot.  Detasheet 
explosive  is  conceptually  more  desirable  than  Primacord  because  it  comes 
in  sheet  for*  and  can  be  rolled  to  any  desired  thickness  greater  than 

15  mils,  the  limit  for  reliable  detonation.  Howeyer.  Detasheet  is 

*  ♦ 

Primacord  is  a  registered  trademark  of  the  Ensign  Bickford  Company. 


(a)  WOVEN  PRIMACORD 


(b)  PARALLEL  STRANDS 


FIG.  8  EXPERIMENTAL  ARRANGEMENT  FOR  FLAT  CHARGE 
FIELD  EXPERIMENTS 
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(c)  l-FOOT  STANDOFF 


FIG.  8  (Concluded) 
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is  relatively  expensive  (about  $10  per  pound  for  Detasheet  C  and  $25  per 
pound  for  Detasheet  D)  and  must  be  ordered  several  months  in  advance. 
Detasheet  was  therefore  used  only  for  equivalent  sheet  thicknesses  and 
standoffs  too  small  to  be  obtained  with  Primacord.  The  very  lightest 
charges  were  made  using  Detasheet  D  extruded  into  0.028-inch-diameter 
rods  weighing  about  4  grains/ft. 

The  majority  of  the  charges  were  made  of  parallel  strands  of 
Primacord  as  in  Fig.  8b  with  a  single  cross  strand  (at  the  left  in  the 
figure)  for  detonation.  Early  in  the  program  a  few  charges  were  made 
by  weaving  Primacord  as  in  Fig.  8a  to  increase  detonation  reliability. 
These  are  time  consuming  to  make  and  the  single  strand  method  of  detona¬ 
tion  proved  just  as  reliable  so  woven  charges  were  discontinued.  In  the 
parallel  strand  configuration,  the  charge  weight  per  unit  area  is 


w  (psf)  = 


12  in. /ft 
7000  grains/ lb 


g(grains/ft? 
S(in. ) 


(18) 


where  g  is  the  Primacord  size  and  S  is  the  strand  spacing.  For  a 
given  charge  density  w,  the  Primacord  size  and  strand  spacing  can 
generally  be  chosen  small  enough  that  the  charge  appears  as  a  continuous 
sheet  at  the  model  standoff.  No  detailed  investigation  of  the  maximum 
allowable  strand  spacing  was  made.  Instead,  a  spacing  rule  was  adopted 
that  allowed  reasonable  shot  assembly  and  then  a  few  shots  were  fired 
at  coarser  spacing  for  comparison. 


The  spacing  mile  is 


where  x  is  the  standoff  distance  between  the  charge  and  the  closest 
point  on  the  model.  This  rule  insures  that  the  strength  of  the  cylin¬ 
drical  shocks  from  individual  strands  at  radius  S/2  where  they  first 
intersect  (see  Fig.  9),  is  large  compared  to  the  overall  "plane"  shock 
arriving  at  standoff  x.  From  Eq.  (43)  for  V  =  2,  the  cylindrical  shock 


pressure  at  S/2 


P 


s 


within  the  self-similar  approximation  is 

8  E2  =  8  .  G 

02(4)2(Y  +  1)  (S/2)2  Y  +  1  4Ci2S 


as) 
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MODEL 


FIG.  9  SHOCK  RADII  FROM  INDIVIDUAL 
STRANDS  AND  FROM  CHARGE 
AS  A  CONTINUOUS  SHEET 


where  G  is  the  energy  release  per  unit  length  of  Primacord  explosive. 
Similarly,  the  plane  shock  pressure  at  x  is 


in  which  we  have  used  =  G/S.  The  pressure  ratio  is 


(21) 


For  x/S  »  16,  p  /p  =  39, 

B  X 

Another  interpretation  of  fixing  x/S  at  a  large  value  is  that 
it  makes  the  pressure  p^  arriving  at  x  from  each  strand,  if  it  were 
to  act  individually,  small  compared  to  the  effective  overall  plane  shock 
Within  the  eelf-siailar  approximation,  p„  is  found  by  replacing  S/2 
in  Eq.  (19)  by  x: 


8 

Y  +  1 


G 

i6v2 
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Then 


9ais  S 

rr~~  =  0.61  - 

lOctgX  x 


(22) 


For  x/S  =  16,  p_/p  =  0.033,  This  result  is  changed  only  slightly  if 
the  more  accurate  expressions  in  Eq.  (60),  taking  into  account  the 
ambient  pressure  p  ,  are  used  to  calculate  p  and  p  .  Many  «experi- 

X  a  X 

ments  were  also  run  with  x/S  =  8,  and  smoothly  decaying  model  pressure 
pulses  were  recorded.  It  was  concluded  that  x/S  =16  is  a  conserva¬ 
tive  spacing  rule  and  that  spacings  at  least  as  coarse  as  x/S  =  8  can 
also  be  used  where  the  minimum  available  explosive  size  does  not  allow 
x/S  =  16. 


Figure  10  gives  a  loading  chart  for  x/S  =  16  and  currently 
available  Primacord  sizes.  The  points  indicate  configurations  for 
which  experiments  were  run  in  the  current  program.  Since  50  and  25 
grain/ft  sizes  have  only  recently  become  available,  load  points  for 
these  sizes  were  run  using  40  grain/ft  Primacord  at  x/S  =  20  and 
x/S  =  10.  Charge  densities  used  ranged  in  multiples  of  2  from  0.0143  psf 
to  0.456  psf.  Corresponding  Detasheet  thicknesses  for  the  same  amount 
of  PETN  can  be  calculated  using  a  PETN  density  of  1.12  gm/cm3,  which  is 

3 

30  percent  of  the  total  1.4  gm/cm  Detasheet  density.  The  resulting 
Detasheet  thicknesses  range  from  2.5  to  78  mils.  It  is  probable  that 
the  coarser  strand  spacing  rule  x/S  =  8  could  also  be  used,  extending 
the  lower  limit  of  available  charge  densities  by  shifting  the  curves 
down  by  a  factor  of  2. 

b.  Charge  Size 

In  the  foregoing  theory,  it  has  been  assumed  that  the  flat 
charge  is  infinite  in  extent.  To  estimate  the  effect  of  finite  charge 
size  on  peak  pressure  and  impulse,  a  series  of  experiments  were  run  using 
rectangular  charges  as  in  Fig.  S  with  sides  of  lengths  b  in  the  ratio 
b/x  =  1,  2,  and  3  relative  to  the  standoff.  Three  series  were  run,  one 
with  x  =  2  feet  and  w  =  0.114  psf,  one  with  x  =  2  feet  and  w  =  0.456 
psf,  and  one  with  x  =  4  feet  and  w  =  0.114  psf.  The  averages  of  the 
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SPACING,  S  -  Inch 


a 


STANDOFF,  x  -  ft 


M-INI-II 


FIG.  10  PRIMACORD  LOADING  CHART  FOR  x/S  •  16  (points 
ore  shots  fired  in  preset  program) 


data  from  all  three  sets  of  experiments  are  plotted  In  Fig.  11,  with 
pressure  and  impulse  (at  0  =  0°,  facing  the  blast)  normalized  to  their 
values  for  b/x  =  2.  The  data  show  that  the  increase  in  pressure  and 
impulse  in  going  from  b/x  =  2  to  b/x  =  3  is  sms  11,  and  suggest  that 
even  at  b/x  =  2  the  pressure  is  within  about  10  percent  of  its  infinite 
sheet  value  and  Impulse  is  within  about  30  percent.  Both  decrease 
rapidly  as  b/x  falls  below  2.  In  the  remainder  of  the  experiments 
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FIG.  II  PRESSURE  AND  IMPULSE  VERSUS 
CHARGE  SIZE/STANDOFF 


♦ 

b/x  was  therefore  fixed  at  2.  Comparison  of  the  absolute  magnitude?  of 
the  pressures  and  impulses  with  the  infinite  sheet  theory  is  deferred  un¬ 
til  later. 

c.  Shock  Path  and  Model  Orientation 

In  addition  to  the  charges  being  of  finite  size,  they  also 
differed  fron  the  simple  one-dimensional  theory  because  they  were  not 
detonated  simultaneously  over  their  entire  area.  Instead,  a  single 
detonator  was  placed  in  one  corner  (see  Fig.  8)  and  detonation  took 
place  along  a  diagonal  line  sweeping  over  the  charge.^  At  some  distance 
from  the  charge  the  shock  angle  is  sufficiently  shallow  that  the  one- 
dimensional  theory  is  adequate.  Figure  12  shows  the  apparent  shock  paths, 


♦ 

At  x  =  1  with  large  charge  densities,  a  diagonally  detonating  charge 
gave  pressures  lower  than  predicted  and  also  a  humped  pressure  profile 
because  of  the  diamond  shape  of  the  charge.  Data  from  these  charges 
were  therefore  deleted  in  favor  of  the  sheet  explosive  charges  to  be 
described  next. 

Detonation  proceeds  at  7  mm/Ps  along  each  Primacord  strand,  parallel 
to  the  side  cf  the  charge.  Since  each  successive  strand  is  initiated 
at  a  slightly  later  time  by  the  edge  strand  running  at  right  angles, 
the  effective  detonation  front  moves  parallel  to  a  diagonal  at  velocity 
7  cos  45  =5  mm/Peec. 
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FIG.  12  SHOCK  TRANSIT  TIMES  AND  APPARENT  SHOCK  PATHS  FOR  DIAGONALLY 
DETONATING  PRIMACORD  CHARGES 

obtained  from  accumulated  data  on  shock  transit  time  to  the  model  from 
the  time  the  detonation  passes  the  center  of  the  charge.  These  are  not 
strictly  the  shock  paths  because  in  each  shot  the  measurement  was  made 
after  the  entire  charge  had  detonated — we  do  not  have  a  steady  detonation, 
so  time  and  distance  cannot  be  interchanged.  However,  since  Fig.  11  indi¬ 
cates  that  the  peak  pressure  is  near  the  infinite  sheet  pressure,  the 
curves  in  Fig.  12  give  a  good  indication  of  the  shock  paths.  Shock  angles 
at  the  model  range  from  about  25  degrees  at  small  distances  from  heavy 
charges  to  about  10  degrees  at  large  distances  from  light  charges. 

Two  orientations  of  the  model  relative  to  the  shock  path  were 
studied,  transverse  and  longitudinal,  as  sketched  in  Fig.  13.  Each 
method  has  advantages.  In  the  transverse  placement  the  shock  arrives 
simultaneously  along  the  length  of  the  model  at  normal  incidence  to  an 
angle  0q,  giving  the  maximum  possible  peak  pressure  and  minimizing  end 
effects.  In  the  longitudinal  orientation  the  shock  arrives  at  angle  9q 
relative  to  the  model  axis,  giving  lower  peak  pressure  and  aggravating 
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(a)  TRANSVERSE  (b)  LONGITUDINAL 


FIG.  13  SHOCK  IMPINGEMENT  FOR  TRANSVERSE  AND  LONGITUDINAL  MODEL 
ORIENTATIONS  {both  drcwn  for  x=  2  feet  and  w=  0.114  psf) 


possible  longitudinal  variation  of  pressure  because  of  end  effects. 

These  disadvantages  were  found  to  be  minor  and  the  advantage  of  knowing 

that  the  pressure  distribution  is  symmetric  about  0  =  0,  'irithout  prior 

knowledge  of  0q,  was  considered  to  be  more  important  in  the  present 

program.  The  majority  of  the  experiments  were  therefore  run  in  the 

longitudinal  orientation  to  avoid  having  to  adjust  the  pressure  gage 

orientation  to  8=0  (the  slopes  0  of  the  curves  in  Fig.  12  were 

o  o 

not  known  at  the  outset  of  the  program).  Longitudinal  detonation  also 
gives  a  continuous  transition  to  the  contact  techniques,  which  also  use 
a  detonation  sweeping  along  the  model  axis.  In  studying  methods  for 
modifying  pressure  distribution,  longitudinal  detonation  is  almost  im¬ 
perative  to  insure  symmetric  pressure  distributions. 

d.  Sheet  Explosive  Charges 

Charges  made  from  Detasheet  were  used  only  at  small  standoffs, 
at  which  the  range  of  equivalent  sheet  thicknesses  available  from  Prima- 
cord  is  limited  (see  Fig.  10).  To  reduce  the  shock  angle  as  much  as 
possible  at  these  standoffs,  Detasheet  charges  were  oriented  with  their 
sides  parallel  to  the  model  axis  as  in  Fig.  14  so  that  the  detonation 


31 


FIG.  14  DETASHEET  CHARGE 


wave  swept  along  the  charge  at  the  full  detonation  velocity  of  7.3  mm/ 
Usee.  Also,  at  small  standoff s, the  size  of  the  model,  rather  than  the 
magnitude  of  the  standoff,  dominates  charge  size  and  b/x  must  be  made 
larger  than  2.  For  0.5  foot  and  1-foot  standoffs,  the  charges  were 
made  with  b/x  =  3. 

Figure  15  is  the  loading  chart  used  for  sheet  explosive.  For 

sheet  thicknesses  greater  than  30  mils,  solid  sheets  of  Detasheet  C  can 

be  used.  At  smaller  thicknesses  Detasheet  D  must  be  used  because  Deta- 

shaet  C  does  not  detonate  well  below  30  mils.  For  the  same  reason, 

below  15  mils  Detasheet  D  was  cut  into  strips  and  spaced  apart  with  a 

spacing  rule  of  x/S  =  12  where  possible,  and  with  x/S  =  6  for  smaller 

average  thicknesses.  For  still  lower  loads,  38-mi 1-diameter  extruded 

rods  were  used;  the  lowest  loads  were  obtained  using  28-mil-diameter  rods, 

the  smallest  cross  section  that  will  detonate.  At  a  standoff  of  3  inches 

* 

the  smallest  obtainable  average  thickness  is  about  one  mil. 


Pressure  and  impulse  data  from  the  3-iuch  standoff  shots  are  not 
reported  because  they  were  not  extensive  enough  to  allow  comprehensive 
treatment . 
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THICKNESS  -  mil* 


FIG.  15  SHEET  EXPLOSIVE  LOADING  CHART  (strip*  ere  oil  15-miJ-thidc 
Datasheet  D;  points  are  shots  fired  in  current  program) 
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To  investigate  methods  for  modifying  the  pressure  distributions 
froa  those  obtained  with  the  flat  charges  described  above,  curved  charges 
and  charges  of  ncmuuiform  density  were  also  used.  These  are  described 
in  Sections  IV  and  V. 

3.  Surface  Loads  froa  Flat  Primacord  Ch 


In  this  section  we  will  examine  the  pulse  shape,  pressure  distri¬ 
bution,  peak  pressures,  and  impulses  froa  Priaacord  charges  at  x  =  1, 
2,  and  4  feet.  In  the  majority  of  the  experiments  the  nodel  cylinder 
was  steel,  6  inches  in  diameter  and  15  inches  long  with  a  1-inch  wall 
thickness.  Pressure  was  measured  with  Kistler  model  601H  and  603H 
gages.  To  neasure  circumferential  pressure  distributions,  gages  were 
mounted  at  the  center  of  the  model  at  22.5  degree  intervals  over  one- 
half  of  the  circumference  and  at  45  degree  intervals  over  the  other 
half.  Possible  axial  variations  in  pressure  were  measured  with  gages 
mounted  3  inches  in  from  each  end,  giving  3  gages  in  all  at  0  =  0 
(facing  the  blast).  These  are  denoted  by  0^,  0^^,  and  C>2,  in  the 
order  in  which  the  blast  sweeps  over  the  model  as  seen  in  Fig.  13b. 

In  Section  III. 6  data  are  given  from  3.5-  and  12-inch-diameter  models 
iastxusented  in  a  similar  manner.  Further  details  of  the  siodels  and 
pressure  instrumentation  are  given  in  Reference  1. 

a.  Pulse  Shape 


As  in  the  pressure  pulses  measured  fro*  explosive  spheres, 
surface  pressure  pulses  fro*  sheet  charges  are  very  nearly  exponential 
in  shape.  Example  pulses  are  given  in  Fig.  16  for  x  =  4  feet  with  w 
ranging  from  0.0275  to  0.228  psf,  and  in  Fig.  17  for  x  =  2  feet  and 
1  foot  at  the  extremes  in  w.  In  Fig.  18  a  pressure  record  typical  of 
those  in  Figs.  16  and  17  is  plotted  on  seai-log  paper  to  show  the  close¬ 
ness  to  an  exponential  decay.  The  measured  pressure  in  this  record  does 
not  differ  froa  an  exponential  until  it  has  decayed  to  about  1/6  of  its 
initial  value. 


(a)  w  =  0.0275  psf,  SHOT  FP-69 
SWEEP  RATE:  200/isac/cr 


fb)  w  =  0.055  psf,  SHOT  FP-68 
SWEEP  RATE:  100/isec/cm 


. C  )  w  -  0.114  psf,  SHOT  FP-70 
SWEEP  RATE.  100/isec/cm 


(d)  w=  0.228 psf,  SHOT  FP-71 

SWEEP  RATE:  UPPER  IOO^sec/cm 
LOWER  200/isec/cm 


FIG.  16  PRESSURE  PULSES  AT  6  =  0  .,  FROM  FLAT  PRIMmCORD  CHARGES 

ir.  i  a 

AT  x  =  4  FEET  (upper  traces  pressure,  filtered  at  120  kc;  lower  traces 
integrals,  unfiltered) 

Because  of  this  exponential  behavior,  semi-log  pressure  plots 
were  used  to  extrapolate  back  to  the  "true"  peak  pressure  P  (3000  psi 
in  the  example)  from  the  highest  pressure  the  gage  was  able  to  record 
within  the  frequency  band  limitation  of  the  gage  and  signal  filter. 
Direct  and  extrapolated  pressures  were  generally  within  15  percent  of 
tach  other.  Impulse  was  measured  by  electronically  integrating  the 
pressure  signal.  These  are  the  lower  traces  in  Figs.  16  and  17. 
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(c)  w  -  0.0275  psf,  SHOT  FP-65  { d  ?  w  *  0.456  psf,  SHOT  FP-74 

2- FOOT  STANDOFF  (SWEEP  RATE  100/isec/cm) 


FIG-  17  PRESSURE  PULSES  /T  4  =  Omid  FROM  FLAT  PRIMACORD  CHARGES 
AT  x  =  1  FOOT  AND  2  FEET 


b.  Pressure  Distribution 

Figure  19  gives  pressure  pulses  measured  at  various  angles 
around  the  front  face  of  tne  model  for  x  =  4  feet  and  w  =  0.114  psf. 

The  pulse  shape  is  about  the  same  at  each  gage,  and  peak  pressure  falls 
off  approximately  as  cos20.  The  records  from  gages  near  0  =  ±90°  have 
a  different  appearance,  largely  because  the  pressure  is  much  lower  and 
the  signal-to-noise  ratio  is  therefore  much  lower.  Pressures  on  the  back 
side  of  the  model  (not  shown)  do  not  rise  suddenly  with  the  shock,  but 
clie  pressures  are  so  low  that  this  difference  in  pul.se  shape  in  unimpcr- 
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(g  )  9  =  i»0°  (UPPER)  AND  270°  (LOWER) 


FIG.  19  PRESSURE  PULSES  ALONG  AND  AROUND  MODEL  (Shot  FP-70,  x=  4  feet, 

w  =  0.114  psf,  signals  filtered  at  120  kc  and  recorded  at  100  psec/cm  except  as  noted) 
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F!G.  20  FRONT  FACE  PEAK  PRESSURES  VERSUS 
CHARGE  DENSITY  FOR  x  «  4  FEET 


to  about  0  =  80°.  Measured  pressures  further  around  the  model  are 
lower  than  given  by  Eq.  (23),  but  in  either  case  are  too  low  to  have 
a  significant  effect  on  structural  response.  Pressure  distributions 
measured  for  x  =  1  foot  and  x  =  2  feet  are  also  represented  by 
Eq.  (23)  within  the  experimental  error  of  about  ±10  percent. 


FIG.  21  MEASURED  PRESSURE  DISTRIBUTIONS  OVER 
FRONT  FACE  OF  MODEL  FOR  x  =  4  FEET 


c.  Peak  Pressure  and  Impulse 

Plots  of  peak  pressures  and  impulses  at  9=0  are  given  in 
Fig.  22.  Although  peak  pressure  depends  directly  on  both  charge  density 
and  standoff,  impulse  over  this  small  range  of  standoffs  appears  to 
depend  almost  entirely  on  charge  density.  One  would  expect  the  impulse 
for  a  fixed  charge  density  to  be  smaller  for  a  larger  standoff,  but  the 
measuremehts  show  that  any  difference  is  hidden  in  the  small  scatter 
in  the  data  of  Fig.  22;  therefore  a  single  impulse  curve  is  drawn  for 
all  three  standoffs. 


IMPULSE — top* 


The  pressures  in  Fig.  22  can  be  compared  with  the  theory  of 

plane  expansion  from  an  intense  explosion  given  in  Section  VI .  To  make 

this  comparison  the  explosive  yield  of  Primacord  is  assumed  to  be 

1570  c&lories/gm.  This  is  at  the  lower  extreiae  given  for  PETN  (from 

7  —  - 

1570  to  2000  eal/gm)  by  Kinney  and  is  chosen  because  the  FETN  here  is 

packed  at  a  relatively  low  density  and,  furthermore,  the  Primacord 

plastic  and  braid  sheath  tend  to  reduce  the  effective  yield  by  absorbing 

more  chemical  energy  than  they  release.  Assuming  this  yield,  the  energy 
2  * 

release  in  ergs/cm  is 


!  =  Tw  P3firigi-gi r_aL . ifi 

1  eff  pS*  lb  _  .2  4‘ 

J  L  L  (30*5 cm)  J  L 


19  x  10 


7  ergs 


4  (“!)■■*•« 

y  cm  f 


x  10  wgff  (psf) 


where  weff  is  the  effective  charge  density.  For  the  charges  here 
*eff  =  2w  because  they  are  on  the  ground  and  expansion  is  confined  to 
one  direction.^  From  Eq.  (7)  the  characteristic  length  is 


E  3  21  *  I"10  V*  =5 

1  cm  4 

s  =  ~  - - —  =  3.16  x  10  w  „  cm 

1  "l  1.013  ,  106  Us"  rft 

cm 


or,  dividing  by  30.5  cm/ft 


(ft)  =  1040  wgff  (psi) 


Using  s^  from  Eq.  (25),  the  pressures  in  Fig.  22  are  xeplotted 
against  .  R  =  x/b1  in  Fig.  23.  Data  from  all  three  standoffs  fall  on  a 


It  is  unfortunate  that  cgs  units  could  not  be  used  throughout,  but  feet, 
Inches,  and  pounds  were  firmly  entrenched  in  the  shot  makeup  and  hard¬ 
ware. 

Several  shots  were  fired  at  x  =  2  feet  from  charges  freely  suspended  in 
air  to  confirm  that  the  effective  charge  weight  is  doubled  for  confined 
charges. 
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FIG.  23  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL 
PEAK  PRESSURES 


single  carve,  demonstrating  that  the  normalization  to  s^  is  appro¬ 
priate.  The  curve  is  therefore  simply  labeled  "running  detonation"  to 
emphasize  that  this  is  the  most  significant  departure  from  a  truly  one¬ 
dimensional  wave.  At  the  low  pressure  end  of  the  data  (overpressure 
ratio  about  20,  or  about  300  psi),  the  measured  pressures  are  only  a 
few  percent  smaller  than  the  theoretical  curve.  Furthermore,  the 
theoretical  and  experimental  curves  have  the  same  slope.  At  higher 
pressures  (above  about  1500  psi)  the  measured  pressures  begin  to  fall 
awgy  from  the  theoretical  curve  until  at  R  =  0.003  the  measured  pressure 
is  less  than  half  the  theoretical  pressure.  This  difference  is  caused 
partly  by  the  neglect  of  real  gas  and  explosive  behavior  in  the  self- 
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similar  theory,  and  pertly  because  the  shock  angle  caused  by  the  running 
detonation  {see  Fig.  13)  increases  with  decreasing  R,  thereby  giving  a 
larger  difference  from  the  normally  reflected  shock  assumed  in  the  theory . 
Pressures  from  shock  tube  experiments  (to  be  described  later),  in  which 
the  shock  was  always  normally  Incident  upon  the  model,  are  in  closer 
agreement  with  the  theory. 


The  data  from  Fig.  22  are  replotted  in  the  peak  pressure- 

impulse  plane  in  Fig.  24.  The  curves  for  c-ach  standoff  are  roughly 

parallel  to  the  constant  time  lines,  giving  t  «#  40  H-sec  at  x  =  1  foot 

o 

to  t  w  160  M’sec  at  x  =  4  feet.  These  curves  are  discussed  more  fully 
at  the  end  of  Section  III. 5. 


4.  Surface  loads  fro  a  Flat  r>°tasheet  Charges 

Pressure  pulses  were  measured  for  sheet  explosive  charges  at 

x  =  0.5  and  1.0  feet.  Typical  records  are  shown  in  Fig.  25  for  x  = 

0.5  foot.  As  at  the  larger  standoffs  using  Primacord,  these  pulses  are 

exponential  in  shape.  Pressure  distributions  differ  from  those  given 

by  Eq.  (23)  because  the  standoff  is  now  so  small  that  it  is  comparable 

to  the  model  diameter.  However,  the  significant  decrease  in  pressure 

caused  by  the  increase  in  shock  travel  distance  with  increasing  9  can 

be  taken  into  account  by  adjusting  at  each  9  the  reflected  pressure 

p  in  Eq.  (24),  using  the  curves  established  for  p  versus  x. 
r  r 

Pressures  calculated  by  this  procedure  agree  with  measured  pressures 
within  the  precision  of  measurement. 


{ b)  h  =  8.8mils 


FIG.  25  PRESSURE  PULSES  AT  6  =  0  FOR  SHEET  EXPLOSIVE  CHARGES 
AT  x  =  0.5  FEET  (sweep  Tate  20  psec/<-m) 


Peak  pressures  and  impulses  are  plotted  against  average  explosive 
thickness  in  Fig.  26.  Again,  to  the  accuracy  of  the  data,  impulse 
depends  only  on  explosive  thickness.  To  compare  the  peak  pressures 
rdth  theory,  is  calculated  using  1570  cal/gm  of  PETN  content 

{about  80  percent  by  weight)  in  the  Detasheet.  This  is  the  same  yield 
used  for  Primacord  and  should  give  comparable  results,  because  the 
20-perc ant  nonexplosive  binder  ingredients  in  Detasheet  reduce  the  yield 
in  the  same  way  as  the  sheath  in  the  Primacord.  The  release 'energy  per 


FIG.  26  PEAK  PRESSURE  AND  IMPULSE  VERSUS  AVERAGE  THICKNESS 
OF  SHEET  EXPLOSIVE  CHARGES 
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unit  area  for  an  effective  thickness  h 


is  then 
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4.19  x  10 


7 


ergo 

cal 


=  7.37  x  1010  h  ..  (cm) 
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3 

where  p  =  1.4  gm/cm  is  the  overall  density  of  the  explosive  and 
0.8  p  is  the  PETN  content.  The  characteristic  length  is 


ergs 


7.37  X  ID10  heff 

S  (cm)  =  — i  - -£2L_ 

1  P1  1.013  xl06fe 


7.26  x  10  h  _  (cm) 
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For  s,  in  inches  and  h  _ 
1  eff 


cm 

in  mils 


s.  (in.)  =  72.6  h  (mils)  (26) 

i  exi 


Using  heff  =  2h,  the  pressure  curvej  in  Fig.  26  give  a  normalized 

curve  very  nearly  coinciding  with  the  normalized  curve  in  Fig.  23  for 

running  detonation  of  Primacord.  Thus  in  Fig.  23  no  distinction  is 

made  between  Primacord  and  Detashe^t.  Similarly,  the  pressure-impulse 

curves  for  Detasheet  are  merely  an  extension  of  the  curves  in  Fig.  24 

for  Primacord,  extending  pulse  times  down  to  t  ta  25  Usee. 

o 

5.  Surface  Loads  from  Flat  Charges  Confined  in  a  Shock  Tube 

At  standoffs  greater  than  a  few  feet,  it  is  impractical  to  maintain 
one-dimensional  flow  by  keeping  the  charge  dimensions  large  compared  •  o 
the  standoff.  A  better  arrangement  is  to  place  the  charge  in  a  shock 
tube  as  in  Fig.  27  so  that  the  walls  of  the  tube  maintain  one-dimensional 
flow.  Some  less  occurs  because  of  friction  along  the  sides  of  the  tube, 
hut  the  comparison  to  theory  given  later  suggests  that  this  effect  is  net 
serious  for  X/D  <  5  and  causes  a  decrease  in  pressure  of  at  most  40  per¬ 
cent  at  x/D  =■  10  (D  i?  tube  diameter).  It  becomes  impractical  to  use 


HEAVY  CHARGE 
v*  « i.26  psf 
400  groint/ft, 
double  grid 


FIG.  27  ARRANGEMENT  OF  FLAT  CHARGES  IN  SHOCK  TUBE 

the  shock  tube  at  standoffs  less  than  a  few  feet  because  detonation  of 
the  charge  must  be  nearly  simu?taneous  to  initiate  one-dimensional  flow, 
A,t  large  standoffs,  multipoint  symmetric  initiation  is  satisfactory,  but 
at  small  standoffs  it  is  more  convenient  to  use  single  point  initiation 
of  a  running  detonation,  which  is  the  scheme  already  described  in  the 
preceding  sections.  Running  detonation  cannot  be  used  In  a  shock  tube 
because  it  results  in  an  oblique  shock  (Fig.  13)  which  would  reverberate 
against  che  tube  walls. 

Charge  arrangement  and  method  of  detonation  in  the  shock  tube  are 
shown  in  Fig.  27.  The  charges  were  made  of  Primacord  grids  and  were 
simultaneously  initiated  at  7  points  by  a  lead-in  spider  of  40  grain/ft 
Primacord.  The  model  was  placed  at  the  open  end  of  the  shock  tube  and 
the  charges  were  placed  at  standoffs  of  1,  2.E,  5,  10,  and  the  total 
20-foot  length  to  the  closed  end  of  the  shock  tube.  At  the  20-foot 
standoff,  expansion  from  the  charge  was  confined  by  the  bottom  of  the 
tube  to  one  direction,  as  in  the  field  charges  on  the  ground.  At 
smaller  standoffs,  however,  the  charges  were  suspended  on  low  density 
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*  3 

!  (about  2  lb/ft  )  styrofoam  pads  so  the  blast  was  free  to  expand  in  both 

I  directions  from  the  charge.  At  these  intermediate  standoffs  the  reflec- 

|  *  ted  shock  from  the  bottom  of  the  tube  arrived  at  the  model  much  later 

\ 

!  than  the  direct  shock  and  was  always  at  a  much  smaller  pressure.  In  the 

I  .  worst  case,  at  the  10-foot  standoff,  the  model  pressure  from  the  re- 

A  X 

i  fleeted  shock  was  less  than  15  percent  of  the  direct  shock  model  pres¬ 

sure  (see  Fig.  28d),  which  can  safely  be  neglected  in  structural  testing. 

j  a.  Pulse  Shape  and  Pressure  Distribution 

t 

Typical  pressure  pulses  at  each  standoff  are  given  in  Fig.  28, 
Pulses  at  various  charge  densities  for  a  10-foot  standoff  are  given  in 
|  Fig.  29.  The  main  body  in  all  the  pulses  is  exponential  in  shape, 

\  followed  by  a  low  pressure  tail  which  decays  at  a  lower  rate.  In  some 

i  cases  the  exponential  is  preceded  by  an  initially  more  rapid  decay.  In 

* 

all  cases,  except  at  the  1-foot  standoff,  the  entire  decay  in  pressure 
\  is  quite  smooth  as  seen  by  the  integrated  pressures  (the  lower  traces). 

At  the  1-foot  standoff  the  method  of  charge  initiation  does  not  give  a 
sufficiently  plane  shock  wave  and  shocks  reflect  from  the  tube  walls. 
Also,  at  one  foot  the  pressure  measured  at  the  center  of  the  model  was 
as  much  as  twice  the  pressure  measured  at  either  end  gage.  However, 
at  x  >  2,5  feet,  pressures  from  all  three  gages  along  the  model  length 
were  generally  within  a  few  percent  of  each  other  and  showed  no  tendency 
for  the  pressure  at  the  center  of  the  model  to  be  larger.  Variation  in 
pressure  around  the  model  was  again  reasonably  well  approximated  by 
Eq.  (23). 

b.  Peak  Pressures  and  Impulses 

Figure  30  gives  peak  pressures  and  impulses  plotted  against 
charge  density.  Each  data  point  is  an  average  from  the  three  0=0 
gages  along  the  shell.  Data  for  x  =  1  foot  are  omitted  for  the  reasons 
just  discussed.  The  pressures  at  x  =  2.5,  5,  and  10  feet  scale  very 
closely  to  the  one-dimensional  blast  theory  in  Fig,  23,  i.e.,  the  curves 
are  merely  shifted  along  the  w-axis  by  factors  of  2.  The  dotted  curve 
for  x  =  10  feet  is  the  curve  for  x  =  5  feet  shifted  by  2  and  falls  very 
close  to  the  solid  curve,  dravn  through  the  10-foot  data.  One  would 
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(a)  x  *  I  !t,  SWEEP  RAT£  =  50/isec/cm  (b)  x  =  2.5 ft,  SWEEP  RATE  *  200/iSec/cm 


(e)  x  =  20 ft,  w#ff  =  0.52psf  x  2, 
SWEEP  RATE  =  (msec /cm 


FIG.  28  SHOCK  TUBE  PRESSURE  PULSES  AT  STANDOFFS  FROM  1  FOOT 
TO  20  FEET  .all  at  w  ^  0.316  psf  except  w  -  0.52  psf  2  in  (e'J 
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FIG.  29  SHOCK  TUBE  PRESSURE  PULSES  AT  10-FOOT  STANDOFF 

expect  the  10-foot  curve  to  fall  below  the  theory  because  of  friction 
losses  along  the  tube  wall.  The  upper  dat-  point  for  x  =  2.5  feet  is 
slightly  below  its  curve,  but  this  is  attrl  uted  to  poor  records  on 
this  shot.  The  pressure  curve  for  x  =  20  ieet  should  fall  on  the 
x  =  10-foot  curve  if  the  one-dimensional  theory  were  still  applicable, 
because  the  charge  is  at  the  bottom  of  the  tube  and  the  effective  charge 
weight  is  doubled.  Instead,  the  curve  lies  a  factor  of  about  1.6  to  the 
right  of  the  solid  10-foot  curve  and  1.7  to  the  right  of  the  dotted 
curve.  This  wide  divergence  from  theory  is  probably  caused  partly  by 
still  further  wall  friction  losses  at  this  large  standoff  and,  perhaps 
more  important,  because  the  charge  was  actually  a  few  inches  above  the 
bottom  of  the  tube,  thereby  spreading  out  the  initial  front  of  the  shock. 

The  magnitudes  of  the  pressures  in  Fig.  30  were  compared  with  the 
intense  explosion  theory  using  Eq.  (25)  and  Fig.  23  as  for  the  field  shots. 
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FIG.  30  PEAK  SHOCK  TUBE  PRESSURES  AND  IMPULSES  VERSUS  CHARGE  DENSITY 

(charge*  fr3e|y  suspended  except  at  x=  20  feet,  charge  against  bottom  of  tube) 


Data  for  x  =  2.5  feet,  5  feet,  and  the  clotted  curve  in  Fig.  30  for 
x  =  10  feet  all  fall  on  a  single  curve  in  Fig.  23.  This  curve  coincides 
with  the  running  detonation  curve  at  the  row  pressure  end,  and  at  higher 
pressures  lies  closer  to  the  theoretical  curve  than  the  running  detona¬ 
tion  curve,  as  one  would  expect,  because  here  we  have  normal  shock  inci¬ 
dence. 

Impulses  for  x  =  2.5  feet  and  x  =  5.0  feet  are  given  by  a 
single  curve  in  Fig.  30  to  the  accuracy-  of  measurement.  The  impulse 
curve  for  x  =  10  feet  lies  about  15  percent  below  the  5-foot  curve, 
and  the  20-foot  curve  (after  it  is  shifted  to  the  right  by  a  factor  of 
2  to  account  for  effective  charge  doubling)  lies  about  25  percent  below 
the  10-foot  curve.  Thus,  as  in  the  field  shots,  impulse  depends  mainly 
on  charge  density,  but  here  the  standoff  is  becoming  large  enough  that 
impulse  per  unit  charge  density  is  decreasing  measurably  with  standoff. 

A  pressure-impulse  plot  of  the  data  in  Fig.  30  appears  in 
Fig.  31  along  with  data  from  the  running  detonation  field  shots.  The 
;c  -  4-foot  curve  for  the  field  shots  is  omitted  for  clarity  and  because 
in  practice  it  is  more  convenient  to  prepare  shock  tube  charges  than 
the  large  8  x  8-foot  charges  required  in  the  field  at  x  =  4  feet.  The 
P,.I  curves  for  the  field  and  shock  tube  data  fall  nicely  together  into 
a  set  cf  nearly  parallel  curves,  but  it  is  apparent  that  at  a  given  stand¬ 
off,  impulse  is  higher  in  the  shock  tube.  Imagining  an  interpolated 
curve  drawn  for  x  =  2.5  feet  for  the  field  shots,  the  corresponding 
x  »  2.5-foot  shock  tube  curve  gives  pressures  from  30  to  50  percent 
larger.  This  difference  is  attributed  to  pressure  decay  in  the  field 
shots  caused  by  the  finite  size  of  the  charges  (see  Fig.  11). 

Also  drawn  in  Fig.  31  are  critical-  curves  for  damage  of  a  weak 

shell  and  a  strong  shell,  both  taken  from  Reference  1.  The  weak  shell 

ia  a  simple  unsupported  6061-T6  aluminum  shell  with  radius-to-thickness 

ratio  a/h  =  100  ami  length-to-diameter  ratio  L/D  =  1.  The  strong 

shell  has  a/h  =  24  and  is  covered  by  a  Micarta*shell  with  a/h  =  12.  The 

sheet  charge  pressure-impulse  curves  span  the  region  between  these  two 

shells  and  between  the  t  =  10  Hsec  and  t  =  1000  Fsec  lines,  as  desired. 

o  o 

— 

Micarta  is  a  registered  tradename  of  Westinghouse  Electric  Corporation. 
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FIG.  31  PRESSURE-IMPULSE  CURVES  FOR  ASYMMETRIC  LOADING  TECHNIQUES 
(6-Snch-diameter  model) 


The  characteristic  times  t  =  I/P  in  Fig.  31  can  be  compared 

to  the  incident  pulse  time  scales  tfJ  calculated  using  Eq.  (6)  from 

the  strong  shock  theory.  For  w  =  1  psf  in  the  shock  tube  at  x  =  5 

feet,  this  equation,  with  s]L  from  Eq.  (25),  gives  tfl  =  458  Psec. 

From  Fig.  31,  t  for  the  point  x  =  5  feet,  w  =  1  psf,  is  200  Psec. 

As  expected,  this  time  is  smaller  than  t^,  The  ratio  being  tQ/t  = 

0.44.  Further  down  the  curve,  at  x  =  5  feet  and  w  =  0,05  psf, 

t  .  —  1750  Psec  and  t  =  400  Psec,  giving  a  ratio  t  /t_,  =  0.23.  On 
ix  o  o  fl 

the  Datasheet  curve  at  x  =  0.5  foot  and  h  =  10  mils  (h  =  20  mils), 

efi 

Eqs.  (8)  and  (26)  give  t  =  36,3  Psec  and  Fig.  31  gives  t  =26  Psec; 

ix  o 
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PEAK  PRESSURE  —  bar* 


I  the  ratio  is  t  /t„  =  0.71.  These  comparisons  show  that  t._  is  a 

,  o  fl  fl 

crude  but  useful  estimate  for  t  ,  with  t  being  about  half  of  t,, 

o  o  fl 

for  high  pressures  and  decreasing  to  about  one-fourth  ox  t  for  '-he 

lower  pressures  and  larger  standoffs  here.  A  better  estimate  for  model 

,  impulse  I  and  time  t  could  be  obtained  by  calculating  the  reflected 

impulse  1  from  a  flat,  rigid  wall  and  taking  t  =  I  /P  .  This  would 
r  r  r  r 

require  extensive  computer  calculations  of  the  type  in  References  8  and 
9,  but  the  usefulness  of  the  results  would  certainly  justify  the  effort. 

All  of  the  curves  in  Fig,  31  extend  easily  to  loads  greater  thar. 

the  strong  shell  damage  curve,  and  all  but  the  x  ~  0.5-foot  curve  can 

be  extended  below  the  weak  shell  curv,?.  At  small  standoffs  the  explosive 

spacing  rule  S/x  >  6  and  minimum  exploci,ve  rod  size  discussed  in 

Section  III. 2  allows  impulses  only  as  small  as  abouv  1000  taps  at  x  = 

0.5  foot.  However,  impulse  as  low  as  100  taps  can  be  obtained  with 
2  5 

gaseous  explosives.  ’ *" 

6.  Effects  of  Model  Size 

^  The  pressure-impulse  curves  in  Fig.  31  are  applicable  only  to  6- inch- 

diameter  models  because  pressure  decay,  and  therefore  impulse,  depends 
upon  model  size.  In  the  limit  of  aa  extremely  large  model  (or  small 
standoff)  the  impulse  is  maximum,  equal  to  the  reflected  impulse  from  a 
flat,  rigid  wall.  As  the  model  becomes  smaller  and  smaller  (or  '•he 
standoff  becomes  larger  and  larger)  quasi-steady  flow  is  established 
around  the  model  in  a  time  short  compa:  ed  to  the  total  blast  -duration. 

The  pressure  quickly  drops  from  the  reflected  pressure  to  a  decaying 
drag  pressure  and  the  impulse  is  determined  by  drag  flow  rather  than 
by  diffraction  flow.  For  very  long  pulses  th‘  diffraction  flow  is  so 
;  short  compared  to  the  drag  flow  that  the  peak  pressure  is  more  appro- 

!  priatcly  associated  with  the  stagnation  pressure.  Then  pressure  and 

i 

impulse  intensity  nc  longer  depend  on  model  size. 

The  6-inch-model  size  was  selected  for  most  of  the  experiments 
described  in  this  report,  because  it  is  a  convenient  size  for  structural 
models.  However,  smaller  and  larger  structural  models  ar^  often  used 
,  and,  to  give  some  indication  of  the  change  in  impulse  with  diameter,  a 
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(o)  2-FOOT  STANDOFF 
w= 0.057  psf 


(b) 4-FOOT  STANDOFF 
w"  0.114  psf 


F)G.  32  ARRANGEMENT  FOR  DETERMINING 

OF  MCDEl  SIZE  (3.5-i«icn  ond  l2-inch-oic<*’>®*er 

cylinder!) 


of  shots  at  x  =  2  feet  are  given  in  Fig.  33  along  with  impulse  data 
from  ihe  previous  shots  on  the  6-inch  model.  The  difference  in  impulse 
on  the  models  is  largest  at  the  smallest  charge  density,  with  impulse  on 
the  32-inch-diametcr  model  being  about  1.4  times  that  on  the  3,5-inch- 
dianeter  model  at  w  =  0.03  psf.  At  the  other  extreme,  w  =  0.45  psf, 
the  impulse  ratio  is  only  about  1.1.  In  a  shot  at  w  =  0.114  psf  and 
x  =  4  feet,  the  impulse  ratio  was  1.3.  These  differences  in  impulse 
are  not  much  lai’ger  than  possible  errors  of  10  to  20  percent  in  reading 
the  impulse,  caused  by  inaccuracies  in  gage  sensitivity  and,  more  im¬ 
portant,  in  low  pressure  drifu  due  to  gage  heating  and  other  effects. 
Nevertheless,  the  smallness  of  the  change  in  Impulse  with  model  size 
observed  here  shows  that  the  pressure-impulse  curves  in  Fig.  31  can  be 
used  as  estimates  for  larger  and  smaller  models ;  more  accurate  estimates 
can  fce  made  over  a  limited  range  by  correcting  with  the  aid  of  the  ratios 
above. 
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FIG.  33  IMPULSES  ON  3.5-,  6-,  AND  12-INCH-DIAMETER  MODELS 
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SECTION  IV 


SIMULATION  OF  SYMMETRIC  LOADS 

The  methods  described  here  for  simulating  symmetric  loads  are 
logical  extensions  of  the  methods  just  described  for  asymmetric  loads. 

To  obtain  quasi-impulsive  loads,  the  flat  charges  on  the  ground  were 
replaced  by  cylindrical  charges  completely  surrounding  the  model,  again 
at  standoffs  of  a  few  feet.  To  obtain  longer  duration  loads,  the  shock 
tube  was  used  as  it  was  for  asymmetric  loads,  but  the  model  was  pointed 
down  the  tube  with  its  axis  ir.  line  with  the  axis  of  the  shock  tube. 

The  experiments  reported  here  are  not  as  extensive  as  those  described 
for  asymmetric  loads  fcfcause  emphasis  was  on  feasibility  of  the  tech¬ 
niques  rather  than  extensive  calibration.  Considerable  further  work, 
both  theoretical  and  experimental,  is  necessary  to  exploit  the  tech¬ 
niques  described  here. 

1.  Surface  loads  from  Cylindrical  Charges 

Cylindrical  charges  2  feet  and  4  feet  in  diameter  were  made  by 
stringing  Primacord  strands  between  cartwheel  end  frames  as  shown  in 
Fig.  34.  The  rigid  model  was  6  inches  in  diameter,  giving  corresponding 
standoff  distances  of  0.75  and  1.75  feet.  For  both  diameters  the  dis¬ 
tance  between  end  frames  was  3  feet  and  the  charges  were  initiated  by  a 
cone  of  Primacord  strands  at  one  end,  resulting  in  a  nearly  circular 
detonation  front  sweeping  axially  along  the  model. 

Pulse  shape  was  generally  exponential,  as  it  was  for  asymmetric 
loads.  Pressure  pulses  at  the  lightest  and  heaviest  charge  densities 
used  at  each  diameter  are  given  in  Fig.  35.  As  Fig.  36  shows,  the 
distribution  of  peak  pressure  around  the  model  was  satisfactorily  uni¬ 
form  in  all  but  one  shot.  In  all  of  the  2-foot-diameter  shots  peak 
pressure  was  uniform  within  ±10  percent.  At  4  feet  the  spread  was 
slightly  larger,  but  still  acceptable,  except  at  the  heaviest  charge 
density,  for  which  peak  pressure  varied  from  4500  to  &000  psi.  Impulse 
was  always  very  nearly  uniform  around  the  model,  even  in  this  shot. 


(a)  2-FOOT  DIAMETER  CHARGE ,  w  =0.0915  psf 
40-grain  PRiMACORD  ON  0.75-inch  CENTERS, 
SHOT  C-88 


(b)  4-FOOT  DIAMETER  CHARGE,  w=  0.023  psf 
40-grain  FRIMACORD  ON  3- inch  CENTERS, 
SHOT  C-S20 


FIG.  34  CYLINDRICAL  CHARGES  FOR  SYMMETRIC 
QUASI -IMPULSIVE  LOADS 


(a)  w  -  0  023  osf  (SHOT  C-42)  (b)  w  =  0.092  psf  (SHOT  C-88) 

2-FOOT  DIAMETER  CHARGES  (x  =  0.75  ft) 


(C)  w  =  0.029 psf  (SHOT  C-86)  (d)  w  =  0.I)5psf  (SHOT  C-87) 

4 -FOOT  DIAMETER  CHARGES  (x  =  1.75  ft  > 


FIG.  35  PRESSURE  PULSES  FROM  CYLINDRICAL  CHARGES  [filtered  afl  20  ^c, 

50  )i sec/cm  sweep  in  (a)  and  (b),  lOOpsec/cm  in  (c)  and  (d)j 

Thus  the  pressure  pulse  at  9000  psi  decayed  more  rapidly  than  the  pulses 
at  4500  psi.  One  would  expect  the  pressure  to  fluctuate  more  for  tne 
larger  diameter  charge  ause  of  imperfect  covergence  and  instability 
of  the  imploding  shock,  out  it  is  not  evident  whether  the  largest  pres¬ 
sure  fluctuations  occurred  at  the  highest  charge  density  by  chance,  or 
because  of  inherently  poorer  convergence  as  charge  density  increases. 
Unsymmet  ic  convergence  also  caused  some  changes  in  the  initial  shape 
of  the  pressure  pulse,  reducing  the  peak  pressure  and  giving  it  a  flat 
or  rounded  top.  These  changes  were  not  serious  and  were  observed  in 
only  a  few  of  the  more  than  50  pressure  records. 
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FIG.  38  PRESSURE-IMPULSE  CURVES  FOR  SYMMETRIC  LOADING  TECHNIQUES 
(6-inch-diometer  model  f) 


on  other  model  sizes  cannot  be  estimated ;  it  would  be  difficult  to  cali¬ 
brate  the  cylindrical  charge  technique  for  a  range  of  model  sizes  with¬ 
out  a  theory  as  a  guide.  Since  the  present  experiments  demonstrate  that 
useful  quasi-impulsive  loads  can  be  obtained  using  cylindrical  charges, 
a  theoretical  program  should  also  be  pursued. 

2.  Symmetric  Loads  in  a  Capped  Shock  Tube 

Longer  duration  loads  were  obtained  using  flat  charges  in  the  2-foot- 
diameter  shock  tube  with  the  model  pointing  down  the  tube  as  shown  in 
Fig.  39..  For  moderately  low  pressures  the  model  can  be  placed  on  a  sting 
with  a  minimum  of  cross  bracing  so  that  the  incident  blast  wave  can  flow 


63 


PEAK  PRESSURE  - bars 


-3R5-  # 


*  -  ^  - - m . . . . *  ,  , . -r ,*ivi&sp«sr&9$5Mfr  rro<$%4?1i£sjE^ldM?'^^  -  ~  ^P«5s-A-^-e 


past  the  model.  Peak  pressure  is  then  near  the  incident  pressure.  At 
high  pressures,  however,  the  necessary  charge  weights  using  the  incident 
pressure  become  objectionably  large.  For  example,  to  obtain  an  incident 
overpressure  ratio  of  200  (5*3000  psi),  Fig.  53  shews  that  the  normalized 
standoff  must  be  less  than  R  =  x/s^  =  0.001.  At  a  standoff  x  =  10  feet, 
Eq.  (25)  gives  a  corresponding  charge  density  of  w  =  x/1040  R  r,  10  psf. 
For  the  2-foot-diameter  shoek  tube  this  gives  a  total  charge  weight  of 
30  pounds.  The  actual  charge  weight  needed  is  probably  about  twice  this 
because  at  these  high  pressures  the  self-similar  theory  overestimates 
the  pressure.  It  is  very  likely  that  a  60-pound  disk  of  explosive  would 
rupture  the  1-inch-thick  steel  wall  of  the  shock  tube.  Sandia  Corporation 


I 

i 


has  obtained  incident  pressures  of  3000  psi  using  plywood  shock  tubes 

10 

which  are  destroyed  with  each  shot. 


To  circumvent  tnese  problems,  symmetric  loads  were  obtained  here 
by  using  the  reflected  pressure  rather  than  the  incident  pressure.  This 
is  done  by  placing  the  model  against  a  rigid  back  plate  as  in  Fig.  39a. 

|  -  The  incident  blast  reflects  from  this  plate  so  that  the  model  is  exposed 

|  to  the  reflected  pressure.  Pressure  records  from  this  configuration  with 

x  =  18.5  feet  are  given  in  Fig.  40.  The  model  is  exposed  to  a  short  part 
of  the  incident  pressure  as  well  as  the  reflected  pressure,  but  the  inci¬ 
dent  pressure  is  much  lower  (about  one-fourth  that  of  the  reflected 
pulses  in  Fig.  40)  and  the.  duration  is  quite  short  (the  time  for  the 
shock  to  traverse  the  model  and  return  from  the  back  plate).  At  the 
pressure  gage  location,  5.2  inches  from  the  back  plate,  this  time  is 
0.5  msec  compared  to  a  total  pulse  duration  of  about  3  msec.  The  impulse 
under  the  incident  pressure  toe  is  negligible,  as  is  seen  in  the  inte¬ 
grated  pressure  traces  in  Fig.  40. 

j  Peak  pressures  and  impulses  from  these  pulses  are  plotted  in 

j  Figs,  41  and  38.  Figure  33  shows  that  a  peak  pressure  of  2000  psi  can 

J  *  be  obtained  with  a  charge  density  of  only  2  psf.  The  lowest  pressure 

j  obtained  here  was  about  300  psi,  but  pressures  as  low  as  50  psi,  the  ex- 

j  *  treme  of  current  interest,  can  easily  be  obtained  with  smaller  charges. 

{  Therefore  no  apparatus  was  built  and  no  tests  were  run  with  the  model 

I 

j  supported  on  a  sting. 


(a ) 


w  =  0.(33  psf  (SHOT  CT-II5 


t 


(b  )  w  ^  0  53 psf  (SHOT  CT— 1 1) 


(c  )  w  =  l.06psf  (SHOT  CT  —  IJ) 


HIG.  40  INCIDENT  AND  REFLECTED 

PRESSURE  PULSES  IN  CAPPED 
SHOCK  TUBE  (sweep  rafe  1  msec,  cm) 
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FIG.  4 1  PRESSURE  AND  IMPULSE  VERSUS  FLAT  CHARGE 
DENSITY  IN  CAPPED  SHOCK  TUBE 


Feak  pressures  are  lower  here  thtn  the  lateral  load  pressures  in 
Fig-  30  because  the  initial  spike  occurs  only  at  the  back  plate  and  is 
lost  by  the  time  the  reflected  shock  has  traversed  back  to  the  gage 
location.  The  impulses  here  art?  about  1.5  times  the  lateral  load  im¬ 
pulses  in  Fig.  30  (again  suggesting  that  calculation  of  reflected  pres¬ 
sure  pulses  will  give  a  good  estimate  of  peak  lateral  loads) .  Figure  38 
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shows  that  the  characteristic  times  of  the  capped  18.5-foot  shock  tube 
piilses  are  near  1.5  msec. 

3.  Recommendations  for  Future  Work 

Figure  33  also  shows  that  there  is  a  wide  gap  in  pulse  times  from 

t  «  0.1  msec  for  4-foot-diameter  cylindrical  charges  to  t  =1.5  msec 
o  o 

for  the  18.5-fort  shock  tube.  The  most  convenient  method  for  obtaining 
intermediate  durations  would  be  to  use  the  shock  tube  at  smaller  stand¬ 
offs  as  was  done  for  asymmetric  loads.  The  minimum  standoff  for  a  9- 
inch-long  model  is  probably  about  5  feet  because  of  difficulties  in  ob¬ 
taining  a  plane  detonation  and  because  at  smaller  standoffs  the  incident 
pressure  pulse  and  shock  decay  along  the  length  of  the  model  beco'no 
significant.  Assuming  that  pressure  in  the  reflected  shock  drops  off  in 
proportion  to  the  total  distance  traveled,  for  a  0.75-foot-long  struc¬ 
tural  model  and  a  5-foot  standoff  the  pressure  will  drop  about  1.3  percent, 
along  the  length  of  the  model.  This  drop  is  well  within  the  overall 
accuracy  of  vulnerability  analyses.  Probably  more  important,  at  small 
standoffs  the  duration  of  the  incident  pressor*  toe  (see  Fig.  40)  will 
become  a  sizable  fraction  of  the  overall  pulse  duration.  Measurements 
of  pressure  pulses  at  various  stations  along  the  model  should  be  made 
at  smaller  standoffs  to  investigate  these  effects. 

Assuming  that  acceptable  pulses  can  be  obtained  at  a  5-foot  stand¬ 
off,  a  line  has  been  drawn  in  Fig.  38  shifted  to  the  left  by  a  factor 
of  4  from  the  18.5-foot  shock  tube  line.  A  sizable  gap  stl.il  remain:' 
between  this  line  and  the  pressure-impulse  curve  for  4-foot~diameter 
cylindrical  charges.  One  possibility  for  providing  pulses  is  this  re¬ 
maining  region  would  be  to  make  larger  cylindrical  charges  with  rigid 
plates  at  each  end  to  reduce  the  amount  of  explosive  needed.  Once  again, 
there  Is  a  need  for  a  theoretical  investigation  of  the  diameters  and 
charge  weights  needed  for  various  diameter  models,  and  of  the  otabiU ty 
of  shock  convergence  for  larger  diameters. 


Shock  tube  length  was  decreased  from  20  to  18.5  fett  because  it  was 
necessary  to  pour  in  a  block  of  concrete  to  seal  the  base  of  the  tube 
against  water  leakage. 
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SECTION  V 

loads  of  other  distributions  and  pulse  shapes 
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As  discussed  in  the  introduction,  the  basic  approach  proposed  for 
determining  structural  response  and  vulnerability  to  transient  surface 
loads  is  to  test  structures  with  pressure  distributions  and  pulse  shapes 
that  simulate  the  principle  features  but  not  all  of  the  detail,  of 
actual  loads.  Indeed, the  actual  loads  are  often  not  known  and  the 
number  of  possible  loads  at  various  target  orientations  is  too  large 
to  contemplate  examining  each  in  detail.  Instead,  response  is  examined 
under  a  basic  set  of  asymmetric  and  symmetric  loads  and  then  under 
loads  perturbed  from  these  basic  loads.  If  significant  differences  in 
response  to  the  perturbed  loads  are  observed,  then  v:e  know  that  response 
to  a  more  extensive  set  of  basic  loa  is  must  be  examined.  If  response 
to  the  perturbed  loads  does  not  differ  significantly  from  the  basic 
load  response,  then  we  are  more  confident  that  the  basic  loads  chosen 
contain  the  significant  features  that  influence  response. 


In  this  section  we  examine  a  few  methods  for  changing  pressure 
distribution  and  pulse  shape.  The  experiment  arc  merely  exploratory 
because  an  extensive  study  of  load  variations  should  be  closely  coupled 
with  the  study  of  structural  response.  This  is  the  subject  of  continuing 
research. 

1.  Changes  in  Load  Distribution 

For  impulsive  loads  from  contact  charges  it  has  been  common  practice 
since  the  inception  of  these  techniques  to  obtain  smoothly  varying  dis¬ 
tributions  of  interest  by  varying  the  explosive  thickness,  or  equivalent 
thickness,  around  the  test  object.  A  logical  question  is  whether  similar 
methods  can  be  used  to  vary  the  distribution  in  quasi-impulsive  loads. 
Since  quasi-impulsive  loads  are  obtained  here  with  sheet  charges  some 

1.  distance  from  the  model,  the  shape,  as  well  as  the  thickness  (density) 

in  the  charge,  can  be  varied.  Here  we  will  describe  a  few  experiments 
to  examine  the  influence  of  charge  shape  and  density  variation  on  load 
;  distribution. 
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Two  sets  of  quasi-impuisive  load  distrioution  experiments  were  per¬ 
formed.  In  the  first,  sheet  charges  of  constant  density  were  applied  at 


-C  , 


•j'c  -  ■  -  - 


i ^^KMar^iaaoRctSssfiBaJaa 


made  at  w  =  0.0143  psf  using  25  graln/ft  Primacord  on  3-inch  centers. 

The  charges  were  all  made  6  feet  square  at  a  standoff  of  1.75  feet  (the 

standoff  of  the  4-foot-diameter  cylindrical  charges).  Four  curvatures 

were  used,  from  flat  to  a  semi-cylindrical  form,  with  the  distance  a 

* 

in  Fig.  42c  constant  between  charges.  The  charges  were  initiated  from 
a  central  detonator  at  the  top  in  Figs.  42a  and  b. 

Pressure  distributions  from  the  constant  density  charges  are  plotted 

in  Fig.  43.  The  curve  for  the  flat  charge  is  of  the  same  family  of 

curves  discussed  in  Section  III  and  is  described  approximately  by  the 
o 

cos  8  distribution  in  Eq.  (23)  (the  pressure  in  Fig.  43  is  about  20 
percent  higher  at  0  =  45°  than  given  by  Eq.  (23)).  As  the  curvature 
is  increased  the  pressures  toward  either  side  of  the  model  are  increased 
as  desired,  but  the  converging  shock  at  0=0°  causes  this  pressure  to 
increase  also.  The  net  effect  is  a  more  gradual  pressure  distribution 
than  for  the  flat  charge  but  still,  even  for  the  2.7-foot  radius  charge, 
the  pressure  at  0  ■-  67.5°  is  300  psi,  only  65  percent  of  the  pressure 
at  0  =  0°.  However,  if  this  distribution  is  compared  with  the  distri¬ 
bution  from  the  flat  charge  that  gives  the  same  peak  pressure  (the 
dashed  curve  in  Fig.  43),  the  flat  charge  distribution  is  much  steeper 
(as  also  seen  by  the  increase  in  the  ratio  of  P^/P^  in  Eq.  (23)). 

The  flat  charge  pressure  at  0  =  67.5°  is  137  psi,  only  31  percent  of 
the  8=0°  pressure,  compared  to  65  percent  for  the  curved  charge. 

From  these  results  it  is  concluded  that  significant  variations  in 
pressure  distribution  can  be  obtained  using  variations  in  curvature  alone. 
Even  broader  variations  can  be  obtained  by  varying  both  charge  density 
and  curvature.  Figure  44  gives  si  diagram  of  the  variable  density  charges 
used  to  demonstrate  this  possibility.  The  density  distribution  is  ap¬ 
proximately  linear  from  zero  density  at  the  center  to  a  maximum  density 
at  the  edge.  Two  shots  were  fired  with  this  distribution,  one  flat  and 
the  other  semicircular,  both  at  a  standoff  of  1.75  feet  as  in  Fig.  42. 
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Pressure  records  from  the  semi-cylindrical  charge  at  constant  density 
0.0143  psf  were  lost  because  of  difficulty  in  estimating  proper  oscil¬ 
loscope  settings. 
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The  resulting  pressure  distributions  are  given  in  Pig.  45.  The  distri¬ 
bution  from  the  flat  charge  is  very  nearly  the  same  as  the  distribution 
from  a  uniform  flat  charge  at  the  same  average,  density,  0.057  psf.  The 
di stribution  from  the  curved  charge  is  markedly  different  from  any  of 

the  other  distributions,  with  a  maximum  pressure  near  0  =  ±80°,  at  a 

c 

value  about  40  percent  higher  than  the  pressure  at  0=0.  In  another 
shot,  using  a  line  charge  (i.e.,  a  purely  cylindrical  blast  wave  ex¬ 
panding  from  the  charge),  which  represents  the  opposite  extreme  in 
charge  density  distribution,  the  pressure  distribution  was  essentially 
the  same  as  from  the  flat  variable  or  flat  uniform  charges. 


FUAT  CHARGE  TO  PRODUCE 
>  SAME  PEAK  PRESSURE 
(SEE  FIG.  21) 
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FIG.  43  PRESSURE  DISTRIBUTIONS  FROM  UNIFORM  C.0143-psf  CHARGES 
AT  VARIOUS  CURVATURES 


AVERAGE  DENSITY 
0.057  psf 


The  similarity  in  distribution  from  line  charges  and  from  uniform 
flat  charges  is  to  be  expected  because  the  shock  radius  at  initial 
contact  (1.75  feet)  is  large  compared  to  the  model  radius  (0.25  foot). 
The  similarity  to  the  distribution  from  the  nonuniform  flat  charge  is 
more  surprising  and  apparently  occurs  because  the  standoff,  although 
smaller  than  the  half- -width  of  the  charge,  is  sufficiently  large  to 
allow  the  stronger  shocks  from  the  outer  region  of  the  charge  to  over¬ 
take  and  coalesce  with  the  weaker  shocks  from  the  center  of  the  charge. 
These  results  demonstrate  that  even  at  standoffs  as  small  as  1.75  feet 
from  a  6-inch-diameter  model,  charge  curvature  is  more  important  than 
charge  density  distribution  in  determining  load  distribution. 

Nevertheless,  if  the  variable  density  charge  is  curved,  the 
variation  of  charge  density  then  has  a  significant  effect  on  pressure 
distribution  because  coalescence  into  a  single,  nearly  plane  shock  no 
longer  occurs.  It  is  apparent  from  Figs.  43  and  45  that  with  appropri¬ 
ate  adjustment  of  curvature  and  density,  pressure  distributions  can  be 

2 

obtained  ranging  from  cos  0  over  one  side  (from  a  uniform  plane  shock) 
to  a  uniform  pressure  on  one  side  and  a  rapid  pressure  decay  over  the 
back  side.  Application  of  these  distributions,  along  with  the  uniform 
distribution,  should  be  sufficient  to  establish  the  sensitivity  of 
structural  damage  to  pressure  distribution  in  quasi-impulsive  loads. 

For  quasi-static  loads,  completely  different  loading  schemes 
must  be  used  to  vary  distribution,  since  the  only  variables  open  in 
the  flat  charge  technique  here  are  charge  density  and  standoff.  These 
affect  mainly  peak  pressure  and  impulse  and  have  already  been  used  as 
the  variables  for  producing  the  basic  set  of  loads.  One  loading  scheme, 
which  conceptually  allows  any  desired  pressure  distribution# is  a  set 
of  vented  chambers  surrounding  the  model,  with  radial  and  axial  vanes 
to  allow  separate  adjustment  of  the  pressure  in  each  chamber.  Develop¬ 
ment  of  this  scheme  would  be  relatively  complicated  and  should  await 
the  results  of  the  impulsive  and  quasi-impulsive  load  distribution  tests 
to  see  how  vigorously  it  should  be  pursued. 
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2.  Changes  in  Pulse  Shape 

Pulse  shapes  were  examined  for  shock  tube  loads  to  compare  the 
pulses  obtained  from  three  methods  of  charge  makeup  in  current  use. 

These  are  (1)  flat  charges;  (2)  uniform  charges  of  finite  length,  giving 
a  reservoir  of  gas  at  a  "constant’’  high  pressure,  as  in  conventional 
shock  tubes;  and  (3)  charges  of  finite  length,  but  having  a  variable 
density  to  give  a  decaying  shock. 

The  effect  of  the  length  of  a  uniform  charge  was  studied  first, 
using  the  charges  shown  in  Pig.  46,  placed  at  the  base  of  the  20- foot 
shock  tube.  Charge  length  varies  from  zero  (flat)  to  6  feet,  with 

total  charge  weight  held  constant.  Pressure  pulses  from  these  charges, 

o  * 

at  0  =  0  on  the  6-inch-diameter  model,  are  shown  in  Fig.  47.  Peak 

pressure  decreases  from  2100  psi  for  flat  charge  to  830  psi  for  the 

6-foot  charge,  but  impulse  remains  essentially  constant  at  100  ktaps. 

For  the  6-foot  charge  the  "steady"  drag  pressure  of  about  160  psi  is 

clearly  seen,  extending  from  t  =  1  to  t  =  4  msec  after  the  initial 

/ 

shock.  Also,  at  t  =  2.1  msec  a  jump  in  pressure  occurs.  For  the 
3-foot  charge  a  "steady”  drag  phase  is  apparent,  but  it  lasts  for  only 
about  1  msec  and  is  indistinct.  The  second  jump  in  pressure  occurs 
earlier,  at  about  t  =  0.8  msec.  For  the  1.5-foot  charge  there  is  no 
steady  flow  and  the  pressure  jump  occurs  at  about  t  =  0.25  msec. 
Finally,  for  the  flat  charge  there  is  neither  steady  flow  nor  a  second 
jump  in  pressure;  the  pressure  merely  decays  smoothly  from  2100  psi, 
with  four  small  pips  probably  caused  by  lateral  reverberations  in  the 
tube. 

These  experiments  demonstrate  that  pressure  pulses  from  finite 
length  charges  differ  considerably  from  the  smooth  pulses  expected  from 
most  weapon  effects.  The  source  of  the  second  pressure  jump  is  not 
known  but  it  obviously  depends  on  charge  length.  For  these  charges, 


The  quantity  a  given  in  Fig.  47  is  called  "charge  ratio"  and  is  a 
measure  of  the  local  (shock  tube)  volume  density  of  the  charge  in 
multiples  of  the  density  of  air  in  the  tube. 
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(a)  L=Of*,  STACKED  FLAI 


(b)  L=  1.5ft,  32  STRANDS  (c)  L=3f»,  16  STRANDS 


FIG.  46  SHOCK  TUBE  CHARGES  OF  VARIOUS  LENGTHS  (all  mode 
of  400  grcin/ft  Primacord,  total  weight  2.74  lb  for  each) 
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FIG.  47  SHOCK  TUBE  PRESSURE  PULSES  F0r.  VARIOUS  LENGTH  CHARGES 
OF  SAME  TOTAL  WEIGHT  (.-sweep  rote  1  msec/cm) 

total  impulse  appears  to  depend  only  on  the  total  amount  of  gas  that 
flows  pant  the  model,  independent  of  how  it  is  generated. 

Variable  density  charges  have  been  proposed  for  producing  decaying 
rather  than'  steady  flow  in  shock  tubes.  Schematic  diagrams  of  three 
such  charges  are  given  in  Fig.  48  along  with  the  corresponding  pressure 
pulses.  In  each  case  the  charge  is  made  most  dense  near  the  model,  at 
a  =  2,  with  succeeding  charge  sections  at  a  =  1  and  £  =  0.5.  In  each 
case  the  pressure  is  roughly  exponential  as  desired,  but  the  pulses  are 
not  smooth,  probably  because  of  the  sudden  changes  in  charge  density. 
This  is  most  apparent  in  Fig.  48c  in  which  the  front  cf  the  charge  was 
only  8  feet  from  the  model,  allowing  little  time  for  discontinuities  to 
smooth  out.  In  no  case  is  the  pulse  as  smooth  as  from  simple  flat 
charges.  It  is  concluded  that  while  variable  density  charges  produce 
useful  loads,  flat  charges  arc  superior  and  should  be  used  if  the  shock 
tube  can  sustain  the  high  local  pressure  near  the  charge. 


(a)  TWO  CHARGE  SECTIONS,  TOTAL  LENGTH  7ft 


3.20  lb  TOTAL 

(c)  THREE  CHARGE  SECTIONS,  TOTAL  LENGTH  i?,f; 


F,G*4?  PRESSURE  PULSES  FROM  VARIABLE  DENSITY  SHOCK  TUBE  CHARGES 
(sweep  rote  I  msec/cm} 


SECTION  VI 

THEORY  OF  BLAST  WAVES  FROM  INTENSE  EXPLOSIONS  IN  A  PERFECT  GAS 

1.  Introduction 

Calculation  of  the  diffracted  pressures  around  the  cylindrical 
models  discussed  in  the  preceding  chapters  is  beyond  the  practical  means 
of  even  the  most  advanced  computer  methods.  Examples  of  the  shock  waves 
and  vortices  around  the  cylinder  have  been  examined  experimentally  by 
Bleakney11  and  give  a  graphic  demonstration  of  the  complexity  involved. 
However,  the  peak  pressure  and  impulse  at  the  leading  edge  of  the  cylin¬ 
der  can  be  estimated  if  the  free-field  blast  wave  is  known.  For  blast 
waves  of  moderate  duration,  in  which  the  flow  around  the  cylinder  is 
dominated  by  the  diffraction  phase,  the  peak  pressure  is  the  reflected 
pressure  of  the  incident,  wave  and  the  impulse  is  a  significant  fraction 
(of  the  order  of  1/2  in  many  of  the  experiments  reported  here)  of  the 
reflected  impulse  from  a  flat,  rigid  wall.  Both  of  these  quantities  are 
relatively  easy  to  calculate.  For  blast  waves  of  long  duration  the 
initial  peak  pressure  is  still  the  reflected  pressure,  bit  it  decays  in 
such  a  short  time  relative  to  the  total  blast  duration  that  for  struc¬ 
tural  response  it  is  of  secondary  interest.  For  these  blast  waves  it 
is  useful  to  assume  that  at  each  instant  the  flow  can  be  treated  as 
quasi-steady  at  the  instantaneous  pressure,  particle  velocity,  ar.d  den¬ 
sity  of  the  incident  wave.  The  peak  pressure  to  which  the  reflected 
pressure  decays  is  then  the  drag  pressure  and  the  impulse  can  be  esti¬ 
mated  from  a  quasi-steady  integration  using  wind  tunnel  data  on  steady 

12 

flow  past  cylinders. 

Perhaps  more  important  than  enabling  one  to  estimate  the  absolute 
magnitude  of  surface  pressure  pulses,  these  simple  calculations  give  a 
means  for  establishing  the  trends  to  be  expected  from  diffraction  pulse 
measurements  such  as  those  given  in  the  preceding  sections.  They  can 
show,  for  example,  the  manner  in  which  peak  pressure  and  impulse  will 
vary  with  charge  size  and  distance-to-charge,  and  the  gain  to  be  expected 
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in  going  from  spherical  to  flat  charges.  Such  calculations  are  given  in 
Section  III.  For  symmetrical  implosion  from  a  cylindrical  sheet  of  ex¬ 
plosive  a  relatively  si-ipie  one  space  variable  theory  can  give  the  com¬ 
plete  pressure  pulse. 

The  present  section  treats  blast  wave  theory  in  an  ideal  gas  with 
these  ends  in  mind.  Section  VI. 2  summarizes  the  existing  self-similar 
solutions  for  explosions  in  spherical,  cylindrical,  and  plane  geometries. 
Section  .VI, 3  gives  an  approximate  theory  for  the  nonself- similar  motion 
of  the  shock  wave,  taking  into  account  the  initial  pressure  in  front  of 
the  shock.  Section  VI. 4  gives  suggestions  for  fruitful  extensions  of 
this  effort. 

2.  Summary  of  tSelf -Similar  Theory 

The  theory  of  blast  waves  from  intense  explosions  in  a  perfect  gas 
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can  be  -found  in  advanced  texts  on  fluid  mechanics,  *  but  it  is  so 
fundamental  to  the  interpretation  of  the  present  blast  load  experiments 
that  we  give  it  here  in  brief  form.  To  aid  in  the  calculation  of  blast 
wave  parameters,  pertinent  numerical  results  for  spherical  and  plane 
blast  waves  are  also  given. 

a.  Characteristic  Length  and  Self-Similarity 

The  explosion  is  visualized  as  a  sudden  release  of  a  finite 

energy  E 3  in  a  small  (spherical)  volume  of  gas,  giving  in  the  limit 

an  infinite  energy  density  at  the  mathematical  point  of  the  explosion. 

As  the  gas  expands  away  from  the  point  of  energy  release  it  drives  an 

intense  spherical  shock  into  the  surrounding  undisturbed  gas.  The 

average  energy  per  unit  volume  within  the  expanding  sphere  decreases, 

eventually  becoming  comparable  to  the  energy  density  in  the  undisturbed 
♦ 

ambient  gas,  initially  at  pressure  and  density  Itor  example, 

* 

Most  of  the  gas  and  energy  are  contained  in  a  thin  layer  immediately 
behind  the  shock,  but  the  average  energy  decrease  with  increasing 
volume  still  gives  a  measure  of  the  local  energy  density  decrease.  Also, 
the  apparent  release  energy  E3  is  related  to  the  chemical  energy 
release  through  a  complex  process  not  treated  here.  Erode®  showed  that 
peak  pressures  from  the  self-similar  theory  are  matched  <in  its  range  of 
validity)  to  numerical  calculations  for  spherical  explosions  from  TOT 
by- taking  £3  equal  to  the  explosive  release  energy.  The  same  result 
was  found  by  Goodnan^  using  compiled  data  from  experimental  pressure 
measurements. 


in  a  perfect  gas  with  specific  heat  ratio  y*  the  ambient  gas  has  an 
internal  energy  (per  unit  volume)  of 


(27) 


The  release  energy  per  unit  volume  o  f  the  blast  sphere  at  radius  r  is 


(28) 


These  energy  densities  are  equal  at  radius 


(29a) 


Similarly,  for  a  cylindrical  explosion  from  a  release  energy  E„  per 

r\  * 

unit  length  along  the  axis,  e^  =  Eg/pr*  and  the  energy  densities  are 
equal  at 


(29b) 


Fcr  a  plane  explosion  from  a  release  energy  Ej  per  unit  area,  eb  = 
Ej/2r  and  the  energy  densities  are  equal  at 


(Y  -  1)E 

r  =  - - - 

e  2px 


(29c) 


Thus,  in  each  geometry  a  characteristic  length  iq 


can  be  defined  by 

(30) 


where  v  =  1,  2,  and  3  for  a  plane,  cylindrical,  and  spherical  explo¬ 
sions,  and  it  is  understood  that  the  release  energy  E^  has  appropriate 
units  in  each  case.  From  the  above  examples  we  see  that  rQ  is  the 
blast  radius  (within  a  multiplicative  constant)  at  which  the  release 
energy  becomes  comparable  to  the  initial  energy  of  the  ambient  gas  en¬ 
veloped  by  the  blast.  For  a  one-pound  pentolite  sphere  in  a  standard 
atmosphere,  rQ  =  10.0  feet. 
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If  we  confine  our  attention  to  the  early  motion  (r  <  <  r  ) 

o 

the  energy  density  within  the  blast  is  so  large  compared  to  the  ambient 
energy  density  p^/(Y  "  D  that  the  ambient  energy  can  be  neglected. 
Stated  another  way,  near  the  explosion  the  shock  pressure  is  so  much 
greater  than  p^  that  p^  can  be  neglected.  Under  these  conditions 
the  ambient  gas  is  described  by  alone  and  there  is  no  characteristic 

length.  Consequently,  with  no  length  as  a  basis  of  comparison,  the  blast 
wave  at  any  instant  is  identical  to  the  blast  wave  at  any  other  instant, 
except  for  scale  changes  as  the  sphere  expands-- the  wave  is  self-similar. 
This  similarity  greatly  simplifies  the  mathematical  treatment  because 
the  flow  can  be  described  in  terms  of  a  single  parameter  §  which  is 
an  appropriate  combination  of  r  and  time  t,  rather  than  in  terms  of 
both  r  and  t  independently.  At  later  times,  as  r  approaches  rQ, 
this  simplification  is  no  longer  appropriate  and  the  pressure  p^  ahead 
of  the  shock  must  be  considered.  This  will  be  done  later.  Also,  the 
self-similar  theory  gives  an  infinite  pressure  as  r  -»  0  and  is  there¬ 
fore  inapplicable  at  small  distances  from  a  real  explosion.  For  small 
r  the  mechanism  of  energy  release  and  the  properties  of  real  gases  must 

be  considered.  These  effects  have  been  treated  by  several  authors,  for 
8  9 

example  by  Brode.  ’ 

In  the  remainder  of  this  subsection  we  consider  that  the  ambient 

gas  can  be  characterized  by  its  density  p  alone  so  that  the  blast  is 

1  * 

self-similar.  In  all,  there  are  five  fundamental  parameters  in  the  self¬ 
similar  problem: 


V  Pj.»  Yf  r»  and  t 


There  are  three  fundamental  units:  force,  distance,  and  time.  The 
fundamental  parameters  can  therefore  be  reduced  to  two  dimensionless 
combinations.  These  are  tak'en  as 


(31) 


* 

Any  other  parameter  can  be  expressed  in  terms  of  these. 


The  specific  heat  ratio  Y  is  a  constant  for  any  given  problem  so  that 
§  is  the  single  independent  variable,  demonstrating  the  self-similar 
nature  of  the  problem.  The  dependent  variables  of  primary  interest  are 
the  pressure  p,  density  p,  and  particle  velocity  v. 

b.  Motion  of  the  Shock  Front 

Before  attempting  to  describe  the  complete  flow  within  the 
blast  sobers,  let  us  examine  the  motion  of  the  shock  front  and  the 
values  of  p,  p,  and  v  just  behind  the  shock.  At  each  (and  every) 
instant  the  pressure  wave  appears  as  shown  in  Fig.  49,  in  which  both 
p  and  r  have  been  normalized  to  their  values  p  and  r  at  the 
shock  front.  Since  §  is  the  only  independent  variable,  each  relative 
position  behind  the  shock  corresponds  to  a  value  of  £  and,  in  particu¬ 
lar,  the  position  of  the  shock  itself  corresponds  to  a  unique  value  §  . 


FIG.  49  PRESSURE  DISTRIBUTION 

IN  SELF-SIMILAR  BLAST  WAVES 


« 
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From  Eq.  (31)  the  position  of  the  shock  is  then 


and  the  velocity  of  the  shock  is 


i 


\H-2 


v  +  2 


(32) 


(33) 


The  variables  at  the  shock  front  can  then  be  calculated  directly  from 
the  strong  shock  relations 


2U 

Y  +  1  ' 


2pjU2 

P2  =  y  +  1 


(34) 


For  example, 


the  pressure  behind  the  shock  is 

vrt-2 


P2  = 


8E  e 
o^o 


(v  +  2)^(y  +  1) 


_1 

V 


(35) 


Thus,  using  only  dimensional  analysis,  all  of  the  quantities  at  the 
shock  can  be  determined  within  an  unknown  constant  which  can  be  deter¬ 
mined  by  a  single  experiment.  The  dimensional  analysis  also  shows  that 
if  all  of  the  flow  variables  are  measured  as  a  function  of  time  at  a 
single  radius  from  any  single  explosion,  these  data  can  be  used  to  calcu¬ 
late  the  complete  flow  for  any  other  combination  of  charge  and  radius,  so 
long  as  the  self-similar  approximatiou  is  valid. 

In  the  next  section  it  will  be  convenient  to  use  nondimensional 

dependent  variables  v7,  p7,  and  p7,  defined  using  Eqs.  (33)  and  (34)  such 

♦ 

that  they  are  unity  at  the  shock  front. 


♦ 

This  normalization  and  the  subsequent  self-similar  solution  follow  that 
given  by  Landau  and  Lifshitz,14  based  on  the  original  work  by  Sedov.13 
The  self-similar  problem  was  also  solved  independently  by  Taylor, and 
by  von  Neumann.17 


V 


I 


I 


f 


I  ** 

I 

I 


4  r  „  / 

(2  +  v)(y  +  1)  t  V 


c.  Flow  Behind  the  Shock 

The  equations  for  gas  flow  behind  the  shock,  are 


(36) 


By 

at 


+ 


v  Sy  1  a£  = 

ar  p  3r 


o 


H  ♦  (|!V,  - 0  <37> 

(**'$($)' 9 

The  first  results  from  Newton's  second  law,  the  second  from  mass  con¬ 
tinuity,  and  the  third  from  assuming  the  flow  is  isentropic  in  a  perfect 
gas.  The  solution  to  these  equations  is  simplified  by  observing  that 
the  first  equation  can  be  replaced  by  an  energy  integral  which  can  be 
found  immediately  for  the  self-similar  problem.  This  follows  from  the 
observation  already  made  that  in  the  self-similar  problem  the  energy  in 
the  undisturbed  gas  can  be  neglected  compared  to  the  energy  Eq  released 
by  the  explosion;  the  total  energy  within  the  blast  therefore  remains 
constant.  Furthermore,  since  we  have  a  similarity  flow,  the  energy  of 
the  gas  inside  any  sphere  of  a  smaller  radius,  which  increases  in  such 
a  way  that  %  =  any  constant  (not  only  §q),  must  remain  constant;  the 

radial  velocity  of  this  sphere,  from  Eq.  (31),  is  v  =  2r/5t  for  v  =  3. 

n 

To  express  this  in  terms  of  the  flow  variables,  consider  an  incremental 
expansion  as  shown  in  Fig.  50.  Particles  that  leave  the  sphere  with 


85 


FIG.  50  INCREMENTAL  EXPANSION  OF  AN  ARBITRARY  SPHERE 
WITHIN  THE  BLAST  WAVE 


2  2 

velocity  v  bring  with  them  an  energy  4ra*  d  v  dt  (w  +  v  /2),  where 
the  enthalpy  w  is  the  internal  energy  plus  the  flow  work  per  unit 
mass: 


w  =  s  +  2  _ _ E _ +  £  = 

p  p(Y  -  1)  p  p(Y  -  1) 


-ELL- 


OS) 


If  the  total  energy  within  the  sphere  is  to  remain  constant,  this  must 

be  equal  to  the  energy  in  the  incremental  volume  swept  out  by  the  ex- 

2  2 

paneling  sphere:  4rrr  p  u  dt  (e  +  v  /2).  Equating  these  energies  gives 


w  + 


:)=  \(«*r) 


(39) 


This  same  equation  results  in  all  three  geometries,  since  the  area 
2 

factor  (4rrrJ  here)  divides  out.  Using  u  =  2r/(v  t  2)t  from  Eq.  (33) 


and  the  dimensionless  variables  in  Eq.  (36),  this  energy  balance  gives 


in  all  three  geometries.  Converting  the  last  two  equations  in  (37)  to 
the  primed  dependent  variables,  and  independent  variables  §  and,  say 
t7  =  t  (in  order  to  properly  differentiate  r/t  in  Eq.  (36))yields 


d  log  § 


d  log  p 
d  log  § 


d 

d  log  5 


(loe  J?) 


(V  +  2)(y  +  1)  -  4v7 
2v7  -  (y  +  1) 


Thus  the  three  partial  differential  equations  (37)  have  been  reduced  to 
an  algebraic  equation  (40)  plus  two  ordinary  differential  equations  (41). 

To  solve  these  equations  it  is  convenient  to  treat  v7  as  the 
independent  variable  rather  than  §  because  derivatives  of  §  and  p7 
can  be  expressed  as  functions  of  v7  alone.  These  can  be  integrated  by 
quadrature.  After  laboii^as  manipulation  this  procedure  gives  explicit 
solutions  for  §  and  p7  in  terms  of  v 7  which,  together  with  p7 
from  Eq.  (40),  completes  the  solution  of  the  problem.  The  resulting 
equations  are  elementary  but  lengthy;  they  are  given  in  Ref.  14,  p.  395, 
for  spherical  geometry  and  in  Ref.  13,  p.  219,  for  all  three  geometries. 

The  parameter  §  (Y?  "is.  determined  such  that  the  total  energy 
in  the  gas  is  equal  to  the  energy  Eq  of  the  explosion.  (The  relation 

between  E  and  the  explosive  yield,  however,  remains  unknown  without 

°  9 

a  study  of  the  explosive  gas  flow.  ) 


-/K^) 


where  A  =  4nr  for  spherical  waves,  A  =  2nr  for  cylindrical  waves, 

13 

and  A  =  2  for  plane  waves.  Sedov,  instead  of  finding  from 

Eq.  (42),  defines  an  artificial  energy  E  =  EQ/a  to  be  used  in  Eq.  (32) 
such  that  §  =  1  at  the  shock  front  and  a  is  found  from  Eq.  (32).  The 


—  (  y4-2  ) 

two  parameters  are  related  by  a  — 


In  this  notation 


used  in  place  of  %.  Figure  5.1  gives  Sedov's  curves  for  a  versus  y 


S7 


FIG.  51  THE  RATIO  E/E  =  a  =  Z~{v+2) 
AS  A  FUNCTION  OF  y 


for  all  three  geometries.  For  air  (y  -  1.4),  a  =  1.075,  1.00,  and  0.85 
for  plane,  cylindrical,  and  spherical  geometries,  respectively. 


d.  Incident  Pressure  Pulses  from  Self-Similar  Solution 


Sedov  also  gives  p/p2»  p/p2*  v/v„ 


as  functions  of 


\  -  bebind  the  shock  in  all  three  geometries.  For  use  in  blast 

load  calculations  we  need  to  find  p,  p,  and  v  as  functions  of  time 
at  a  fixed  radius-  rf.  Here  we  find  the  incident  pressure  pulse  p(t), 
the  quantity  of  most  interest  for  estimating  diffracted  pulses  from 
available  data.  From  Eqs.  (35)  and  (32)  the  peak  pressure  is 


8E 


Pf  = 


a(v  +2)  (y  +  1)  r 


v 


(43) 


at  time 


*f  = 


/  \l/2 

(?) 


yf2 

2 


(44) 


"4 
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At  a  later  time  t  the  shock  front  r^  has  continued  past  r  =  r^  and 
the  pressure  at  r^  has  decayed  because  both  p^/p^  (see  Fig-  49)  and 
p  decrease  with  time.  We  define  a  normalized  pulse  time  t  after  the 

a 

shock  has  passed  according  to 


V*-2 


(45) 


in  which  Eq.  (32)  has  again  been  used.  The  pressure  at  t  is 

P(t)  =  p. 


fe)  =  p'(^)  "fe) 


(46) 


Equations  (45)  and  (46)  are  parametric  equations  for  p  and  t  with 
r./r  as  a  parameter.  We  can  now  replace  r  by  r  without  confusion 

I  A  X 

so  that  these  equations  become 


yt-2 

7  =  X  2  -  1  ,  =  XVP  (X)  (47) 

Pf  V 

13 

where  P  (X)  is  p/p  vs.  \  =  r/r  as  tabulated  by  Sedov,  pp.  222- 
V  «  « 

223,  and  given  here  in  Fig.  49. 


Curves  of  p(t)/p.c  are  given  in  Fig.  52  for  v  =  1,  2,  and  3. 

The  general  shapes  of  the  curves  do  not  differ  appreciably  for  small  j, 
but  for  large  t  the  pressure  decays  more  slowly  in  cylindrical  and 
plane  geometry  than  in  spherical  geometry.  The  difference  is  particularly 
evident  for  plane  geometry.  (See  also  Fig.  6  for  a  comparison  with  the 
time  scale  for  the  plane  wave  compressed  to  give  close  agreement  at  early 
times.)  For  t  >  >  1  the  limiting  forms  of  Eqs.  (47),  along  with  the 
observation  that  P^  (X  =  0)  are  constants,  yield 


-♦  P  (0)  T 


2y 

vf2 


(48) 


—2/3  -1 

Thus,  for  large  times  the  incident  pressure  decays  as  t  ,  t  ,  and 

/c 

t  in  plane,  cylindrical,  and  spherical  geometries. 
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(a)  LINEAR  SCALES 
y  *1.4 


FIG.  52  INCIDENT  PRESSURE  PULSE  AT  r  =  rf  FOR 
NORMALIZED  10  p,  =  8E0/a(i/ +  2)2(y  +  1)^,  tf 
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3 .  Shock  Front  Motion  into  an  Initial  Pressure 

la  the  remainder  of  this  section  we  will  treat  blase  waves  taking 
into  account  the  initial  pressure  (and  energy)  of  the  ambient  gas  ahead 
of  the  shock  front.  As  in  the  self-similar  theory,  determination  of  the 
motion  of  the  shock  front  is  relatively  simple  and  will  be  undertaken 
before  the  complete  flow  behind  the  shock  is  studied.  Since  we  include 
now  the  energy  entering  the  blc.st  sphere  as  it  envelopes  the  surrounding 
gas,  the  shock  will  decay  less  rapidly  than  in  the  self-similar  theory, 
in  which  the  total  energy  was  taken  to  be  constant. 

A  method  for  finding  a  first  (and  second)  approximation  to  the 

18 

modified  shock  froat  motion  was  suggested  by  Korobeinikov.  Ke  com¬ 
pared  the  dependence  on  r  of  the  particle  velocity  v  just  behind 
the  shock  front  as  given  by  the  self-similar  theory  (valid  for  small 
r2)  and  by  the  asymptotic  theory  of  Landau  (pp.  375-377,  valid  for 
large  r 2>  and  observed  that  there  was  little  change  in  this  dependence. 
This  is  demonstrated  in  Table  I.  On  the  basis  of  this  observation,  he 
suggested  that  for  intermediate  values  of  r  the  particle  velocity  can 

A 

be  taken  from  the  self-similar  theory  and  then  the  shock  velocity  U  can 
be  calculated  from  the  particle  velocity  jump  condition  across  the  shock. 
Using  this  modified  expression  for  U  the  pressure  p  and  density 

A 

can  then  be  determined  from  the  other  two  jump  conditions. 

Table  I 


DEPENDENCE  OF  PARTICLE  VELOCITY  v  ON  SHOCK  RADIUS  r_ 


Theory 

Geometry 

v  =  1 

V  =  2 

v  =  3 

plane 

cylindrical 

spherical 

Self-similar 

-1 

r2 

-3/2 

r2 

Acoustic  shock 

-3/4 

r2 

-1  -1/2  /  r2  \ 
r’  lcz  [~J 

SI 


Combining  Eq.  (53)  with  the  first  of  Eqs.  (49)  gives  the  desired  ex¬ 
pression  for  shock  velocity: 


M  =  A 


♦  VPT 


(54) 


where 


a2  = 


R 


-v 


ay(v  +  2) 


(55) 


iiut3titt'ting  Eq.  (54)  into  the  second  of  Eqs,  (49)  gives  the  decay  of  the 
shock  overpressure  with  distance. 


-  1  -  -fit 


1  +  (1  +  A"2) 


1/2  ' 
•i 


(56) 


The  limiting  form  for  small  radius  R  (lar^e  A)  is 


l2.  -  1  =  _5l  .2 


v  +  1 


A  = 


8R 


<X(V  +2)  (y  +  1) 


(57) 


After  multiplying  by  p^  and  taking  p^  -»  0,  this  is  seen  to  coincide 
with  the  seif -similar  formula  (35).  The  density  p2  can  be  similarly 
found . 

For  large  R  the  self-similar  particle  velocity  (51)  becomes  a 
poor  approximation  for  cylindrical  and  spherical  waves  (see  Table  I). 

In  this  range  the  small  disturbance  expressions  " able  I  can  be  used, 
matching  the  large  and  small  R  expressions  at  some  intermediate  value, 
For  example,  in  the  spherical  case  Korobeinikov  matches  them  at  R  =  2 
to  obtain 


1/2 

5(y  +  D(ay) 

4R~1[log(R/2)  +  i]"1/2 
5(y  +  i)<2aY)1/2 


R  £  2 


(58) 


P  2  2 


Similarly,  for  cylindrical  blast  waves  we  take 


R  £  2 


(Y  +  lXccY) 


21/4(y  +  IXoly)1^ 


R  a  2 


Using  these  exnressions,  the  incident  overpressure  throughout  the  entire 
range  of  R  is  given  by  Eq.  (56)  with  the  following  expressions  for  A  2 


V  =  1  Plane  waves 
(a  =  1.075  for  y  =  1.4) 


A-2  =  OaYR 


all  R 


v  =  2  Cylindrical  waves 
(a  =  1.00  for  y  =  1.4) 


lOaYir 


i6yrOTR3/2 


R  £  2 


R  £  2 


V  =  3  Spherical  waves 
(a  =  0.P5  for  y  =  1.4) 


25aYR‘ 


A"2  » 


R  5  2 


50aYR2[£n(R/2)  +  l]  R  &  2 


The  reflected  overpressure  ratio  is  given  by 


-  1  = 


(2H-  +  DPg/Pj  +  1 

lVpi + 1 


i-i 

pi 


where  H  =  (y  -  l)/(y  +  1). 

Figures  53  and  54  give  curves  of  incident  and  reflected  overpressure 
ratios  from  these  equations.  Also  shown  are  (dashed)  curves  from  com¬ 
piled  data  on  spherical  pentolite  explosions,  taken  from  Goodman.15  The 
approximate  theory  matches  identically  with  the  experimental  data  for 
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FIG.  53  INCIDENT  SHOCK  OVERPRESSURE  RATIO  vs.  R  =  r/fEo/p,)17" 

R  >  0.2  if  Eq  is  calculated  on  the  basis  of  1430  cal/gin  (pressures 
were  measured  only  out  to  R  =  7).  At  smaller  radii  real  gas  and  ex¬ 
plosive  properties  become  important  and  the  measured  pressures  are 
smaller  than  predicted  by  the  self-similar  theory.  At  larger  radii 
other  mechanisms  of  wave  damping,  such  as  viscosity,  become  important. 
However,  the  range  of  validity  of  the  approximate  theory  is  quite  broad, 
giving  reliable  reflected  overpressure  ratios  from  about  0.1  to  1(>C  a 
useful  range  for  quasi-impulsive  structural  response.  This  excellent 
agreement  for  spherical  explosions  gives  some  confidence  in  using  the 
approximate  solutions  for  cylindrical  and  plane  blast  waves,  for  which 
few  experimental  data  are  available. 


FIG.  54  REFLECTED  SHOCK  OVERPRESSURE  RATIO  vs.  R  =  r2/(Eo/p1)1/*' 

4.  Future  Theoretical  Program 

In  order  to  determine  the  pressure-time  history  behind  the  shock, 

the  complete  flow  must  be  calculated.  Although  extensive  calculations 

8  9  13 

have  been  made  for  spherical  blast  waves,  ’  *  the  authors  know  of  no 

published  results  for  plane  waves  other  than  the  self-similar  curves  of 
13 

Sedov  which  were  used  here  in  Section  VI. 2.  The  usefulness  of  plane 
blast  waves  in  structural  testing,  as  demonstrated  in  the  preceding 
sections,  gives  sufficient  motivation  for  making  plane  wave  calculations. 
The  calculations  should  include  a  complete  description  of  both  the  flow 
behind  the  shock  as  it  moves  out  from  the  explosion  and  the  reflected 
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wave  from  a  flat,  rigid  wall  at  a  sequence  of  values  of  R.  The  reflected 
pressure-time  histories  would  give  a  first  estimate  of  the  diffracted 
pressure-time  history  on  structural  models  and  an  upper  bound  on  the  im¬ 
pulse. 

These  calculations  could  be  made  using  the  self-similar  solution  as 
initial  conditions  at  a  radius  small  enough  that  the  ambient  pressure 

g 

can  be  neglected,  as  was  done  in  the  spherical  problem.  Similarly,  the 
self-similar  theory  could  be  used  to  start  the  calculations  for  the  im¬ 
plosion-explosion  problem  from  cylindrical  sheet  charges.  In  this  case 
the  reflected  pressure  pulse  at  the  rigid  cylindrical  boundary  (model) 
would  be  the  true  model  pressure  within  the  accuracy  of  the  theory.  Such 
calculations  would  allow  relatively  inexpensive  determination  of  peak 
pressure  and  impulse  as  functions  of  charge  diameter  and  thickness  for 
various  model  sizes.  They  should  be  accompanied  by  an  investigation  of 
the  stability  of  the  implosion  process,  which  defines  the  allowable  pa¬ 
rameters  for  acceptably  small  pressure  fluctuations. 
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